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This brief review discusses electronic properties of mesoscopic graphene-based structures. These 
allow controlling the confinement and transport of charge and spin; thus, they are of interest not 
only for fundamental research, but also for applications. The graphene-related topics covered here 
are: edges, nanoribbons, quantum dots, pn-junctions, pnp-structures, and quantum barriers and 
waveguides. This review is partly intended as a short introduction to graphene mesoscopics. 
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I. INTRODUCTION 

Graphene is a two-dimensional (2D) layer of carbon atoms ordered into a honeycomb lattice as shown in Fig. [T] 
It is a material with a host of unusual properties flj-fioj including (among others): Dirac spectrum of low- lying 
quasiparticles Q , large mean- free- path Q , and high electron mobility [U EH . 

Besides its purely fundamental importance, researchers view graphene as a promising new material for electronic 
(ljjj . chemical [l4[, or electromechanical pl| applications, where graphene's unique properties may be of substantial 
benefit. Unlike 3D matter, whose bulk is hidden from direct observation and influence, graphene's "bulk", its 2D 
surface, is always exposed, and its structure may be inspected or modified with greater ease. Furthermore, the Dirac 
energy dispersion in 2D implies that graphene is a gapless semiconductor, whose density of states vanishes linearly 
when approaching the Fermi energy. As such, it is "a bridge material" separating the worlds of semiconductors (with 
an energy gap between the valence and conducting bands) and metals, with a finite density of electronic states at 
the Fermi energy. Depending on the operating regime, graphene can be pushed in either direction. For exa mple, i t 
is possible to open a gap in a sample with the help of chemical modifications 0, [I?} , or lateral confinement (l8l - [2C| . 
Alternatively, one can make graphene metallic, e.g., by chemical doping (2lj . Some graphene samples have spatially- 
varying electronic properties, due to local modifications on the sample. The long electronic mean-free-path, which can 
be of the order of micrometer, implies that electronic signals can travel unimpeded large distances through a device. 
These features might be very useful in applications. 

The unusual properties of graphene motivated significant research efforts. The field grows very fast: the ISI web 
site reports that by October 2010 there were more than 5,000 publications with the word "graphene" in their titles. 
Clearly, this is an enormous volume of scientific work, of which our brief review covers only a very small fraction. Its 
scope is very limited in several respects. As it is obvious from the title, we direct our attention to mesoscopic graphene 
systems, a topic at the boundary between fundamental and applied research. Furthermore, we mainly discuss the 
electronic aspects of graphene mesoscopic systems, especially those which may be relevant for possible electronic or 
spintronic applications, for example, charge/spin transport and confinement, and control over them. Lattice properties 
are dealt with only when the lattice affects the electrons significantly. Several topics are deliberately omitted due 
to space constraints; these include: quantum Hall effect, thermal transport phenomena, phonons, and mechanical 
properties of graphene. 

The review is organized as follows. In Sec. [II] we discuss the most basic electron properties of an infinite graphene 
sheet. The physics of graphene edges is reviewed in Sec. IIII1 Sections IIV[ fVl and IVII focus on nanoribbons, quantum 
dots, as well as pn-junctions and pnp-structurcs, respectively. Sec. IVIIl discusses the barriers created by the combined 
application of magnetic and/or electric fields. Conclusions arc presented in Sec. IVIII1 The main part of the review is 
kept non-technical for it to be accessible by a general reader. More involved discussions are relegated to Appendices. 

II. BASIC PHYSICS OF A GRAPHENE SHEET 

For completeness, in this section we quickly remind the reader the basic single-electron properties of a graphene 
sheet. A more detailed presentation can be found in Appendix [X] It is common to describe a graphene sample in 
terms of a tight-binding model on the honeycomb lattice. Lattice parameters for graphene, as well as some other 
microscopic characteristics, are given in Table |TJ Honeycomb lattice can be split into two sublattices, denoted by A 
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armchair edge 



FIG. 1: (Color online) Geometry of the graphene lattice showing: primitive lattice vectors ai,2, diatomic lattice unit cell 
(dotted-line rhombus), and vectors £1,2,3, connecting the nearest neighbors. Red (black) circles correspond to the A (B) 
sublattice. Three different types of edges (zigzag, armchair, and Klein edge) are shown. The Klein and zigzag edge violate the 
symmetry between the sublattices (the atoms at the edge sites belong exclusively to sublattice A: they are all red), while the 
armchair edge does not (it has both black and red atoms). 



Graphene parameters 


Value 


C-C bond length, ao 


1.4 A 


Lattice constant 


2.46 A 


Hopping amplitudes: 
nearest neighbor, t 
next-nearest, t' 
third-nearest, t" 


2.8 eV 
0.1 eV 
0.07 eV 


Fermi velocity, i>f 


1.1 x 10 b m/s 



TABLE I: Graphene parameters at a glance. 



and B. The Hamiltonian of an electron hopping on a graphene sheet is given by 

H = -tJ2 E (1) 

R.G.A t=l,2,3 

where R runs over sublattice A, and t = 2.8 eV is the nearest-neighbor hopping amplitude. The vectors Si (i = 1, 2, 3) 
connect the nearest neighbors (see Fig. [T] showing the geometry of the graphene lattice). When necessary, H can 
be augmented by interaction or longer-range hopping terms (see Table U for values of the longer-range hopping 
amplitudes) . 

Since there are two atoms in graphene's unit cell, it is convenient to describe the single-electron wave function of 
graphene as a two-component spinor ^ . This introduces an isospin quantum number. For every momentum k lying 
within the Brillouin zone, Fig. [21 the Hamiltonian H has two eigenvalues £k± , which have the same magnitude and 
opposite signs. The eigenvalue £k+ > (ek- < 0) corresponds to the conduction (valence) band of graphene. 

The functions £k± vanish at the six corners of the Brillouin zone: Ki,2 = (0, ±47r/(3-\/3ao)) and K3 4,5,6 = 
(±27r/(3ao), ±27r/(3v / 3ao))- Here the symbol ao denotes the carbon-carbon bond length. Near Kj, i = 1,...,6, 
the dispersion surface can be approximated by two cones with a common apex 

Ekiocilk-Kil. (2) 

The conduction and valence bands touch each other at the cones' apex. 

Of the six cones only two can be chosen to be independent: the remaining four arc connected to these two by a 
reciprocal lattice vector. Thus, the cones Ki,...,6 can be split into two equivalence classes. These classes are commonly 
denoted by K and K', and referred to as 'valleys'. 
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K 

FIG. 2: (Color online) The Brillouin zone of graphene is a perfect hexagon. The Dirac cones are located at the corners of the 
Brillouin zone. The six cones can be split into two equivalence classes (cones within the same class are connected by dashed 
lines). These classes are commonly referred to as K and K' . 



When graphene is not doped, its Fermi level passes through the cone apexes. In such a situation, if one is interested 
in the low-energy description, only the states near the cones must be accounted. For states with energies near the 
cone apexes, it is possible to use the following Weyl-Dirac equations 

Eflflfl = i?*l,2, (3) 

H = -ihv F (a y d x ± <T x d v ) = hv F (g^^. ^ 9x q %Qy ^ . (4) 

which have dispersion as in Eq. ([2]). These equations become invalid away from the cones. The spinor wave function 
^i (^2) corresponds to the electron states near the cone K (K'). The plus (minus) sign in Eq. © corresponds to K 
{K'). The low-energy physics of electrons in graphene is equivalent to four species of two-dimensional massless Dirac 
electrons: two different spin directions and two cones, K and K' , giving overall fourfold degeneracy. 

Pristine undoped graphene is a gapless semiconductor. This means that its density of states does not have a gap, 
but vanishes linearly when the energy approaches the apexes. Sometimes it is desirable to open a gap in the graphene 
spectrum. As shown in Table [TT1 this can be achieved by employing various mechanical, electronic, and/or chemical 
methods. In particular, in monolayer graphene the gap can be induced by substrate or strain engineering [22|, I23L a s 
well as by deposition or adsorption of molecules on the graphene layer, such as, for example, water and ammonia [24J. 
Based on numerical studies, the value of the energy gap can range from a few meV to hundreds of meV. Most 
importantly it can be larger than room temperature as required for graphene-based transistors. In bilayer graphene 
the gap can be induced and continuously tuned, for instance, chemically through selective doping [25l |. or even 
electrically by applying gate voltages [26|. The fact that graphene's band structure can be controlled externally and 
with rather simple processes is a nontrivial result which reveals the potential of graphene for nanotechnology. 



III. EDGES OF GRAPHENE SAMPLES 



The characteristics of a mesoscopic device depend substantially on its edges. Therefore, it is important to study 
the electron behavior near the graphene edge. 
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Inducing an energy gap in graphene 


Method 


Gap in monolayer (meV) 


Gap in bilayer (meV) 


Nanoribbons* (width ~ 15 nm) 


200 




BN-ft / Cu(lll) substrate [27] 


53 / 11 




SiC substrate* [281 


260 




External square superlattice 


65 




Strain engineering [22] 


300 




Adsorption of molecules [30J 


2xl0 :! 




Graphene covered by H 2 / NH 3 [24] 


18 / 11 


30 / 42 


Nanoribbons* (width ~ 30 nm) [311 




50 


Electrical gates* [261 




250 


Selective doping* (potassium) [25] 




100 


Electric field effect* [32] 




150 



TABLE II: Brief summary of possible methods to induce an energy gap in monolayer and bilayer graphene. Asterisks (*) 
indicate experimental demonstrations; otherwise the value of the gap is a theoretical prediction. In some cases, the energy gap 
is tunable and its exact value critically depends on the details of the specific method. Here, 'BN-/i' denotes boron nitride in 
the hexagonal configuration. 




FIG. 3: Reczag, or ZZ 57, edge of a graphene sheet, from [3E|. In the right panel the edge unit cell is shown. It consists of 
a pentagon and a heptagon. The latter polygons are the reason why this edge type is called '57'. Numbers in the right panel 
are the bond lengths in A. Reprinted figure with permission from P. Koskinen, S. Malola, and H. Hakkinen, Phys. Rev. Lett. 
101, 115502 (2008). Copyright (2008) by the American Physical Society. 



A. Edge-stability issues 

Two kinds of edges are often discussed in the literature: zigzag and armchair. A form of the zigzag is the Klein 
edge [33j . All three types are shown in Fig . [T] They are the most symmetric variants of edges in graphene. More 
complicated edges were also studied [34H38j . 

Of course, in a laboratory sample some of these edge types may be unstable chemically or undergo reconstruction. 
The possibility of the edge reconstruction has been addressed in several publications. Most importantly, it appears 
that the pristine zigzag edge is unstable: recently, it was predicted on the basis of density-functional theory (DFT) 
calculations [3!| that it might undergo reconstruction at room temperature, and become a reczag (short for 'recon- 
structed zigzag' [i^, see Fig. [3]). This kind of edge is often called 'ZZ 57'; namely, it is a version of zigzag (thus 
ZZ) edge, in which the edge hexagons are replaced by pentagons and heptagons (hence the 5 and 7). Experimental 
data supporting the existence of the reczag edge were presented in |40j ] . A similar conclusion was reached in [38| : the 
energy of the zigzag edge is substantially higher than the energy of the reczag. In Refs. [4l|, [42j the non-hydrogenated 
zigzag edge was not listed among stable configurations. 



The conclusions of Ref. [39 
where the experiment of Ref. 



regarding the relative stability of zigzag and reczag, were challenged in Ref. [43j, 
L 44j |. proving the existence of the zigzag edge in a laboratory sample, was quoted. The 
experimental demonstration of the zigzag edge stability was also reported in Ref. . Results of Ref. [46| are also 
in disagreement with Ref. [39j]. However, the authors of Ref. (46j were unsure whether their molecular dynamics 
simulations can provide a reliable answer to the question of the edge stability. 

The chemical stability was also investigated. It was pointed out in Refs. [4l|, |42j on the basis of DFT calculations 
that the reczag and armchair are stable only when the concentration of hydrogen in the surrounding media is very 
small. If this is not the case, other types of edges, with hydrogen atoms attached, are stabilized (see Fig. 2]). The 
results of DFT are consistent with Clar's theory of the aromatic sextet [47l l48j. 

The DFT calculations of Ref. [38| demonstrated that the energy of zigzag, armchair, reczag, and more complicated 
regular edges always decreases upon monohydrogenation. This agrees with Ref. j4lj : when enough hydrogen is present 
in the surrounding media, hydrogenation of the edge occurs. 

These results suggest that the edge stability is a complicated problem in graphene. The edge stability depends on 
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k(27t/L) k(27u/L) k (2jt/L) k (2n/L) k (2rc/L) 

FIG. 4: Stable types of graphene edges, from Ref. [4l| . Row A of the figure shows the five most stable configurations of 
the graphene edge with and without hydrogen attached to the unsaturated bonds. There, hydrogen is represented by small 
black circles. Rows B and C show other stable armchair and zigzag terminations. The DFT calculations reported in Ref. [4ll | 
predict that monohydrogenated zigzag denoted as zi in the figure and reczag denoted as z(57)oo are stable only at extremely 
low hydrogen concentrations. At standard atmospheric conditions, an, Z211, and a22 are the most stable types of edges. Note 
that the pristine zigzag edge (studied in numerous papers) is not listed as a stable configuration. More complicated types of 
the graphene terminations are presented in panel D. In panel E the representation of the benzenoid aromatic carbon ring as 
a superposition of two degenerate Kekule configurations is shown. Reprinted figure with permission from T. Wassmann, A. P. 
Seitsonen, A. M. Saitta, M. Lazzeri, and F. Mauri, Phys. Rev. Lett. 101, 096402 (2008). Copyright (2008) by the American 
Physical Society. 



Graphene edges 




Zigzag 


Armchair 


Klein 


Reczag 


Stability [39] 


Unstable 
(to reczag) 


Stable 




Stable 


Edge states 


Yes 


No 


Yes 




Magnetism [411 


Ferromagnetic 


No 




No 


Stress [49] 


Compression 


Compression 




Tension (weak) 



TABLE III: Different properties of graphene edges. In addition to the three types presented in Fig. [JJ (zigzag, armchair, and 
Klein edges), a reconstructed zigzag (reczag) edge is now included in this comparison. 



the orientation of the edge and is affected by the chemical environment. 



B. Electrons near edges 

The simplest way to describe an electron near the edge is to resort to the Wcyl-Dirac equation (|3|) with appropriate 
boundary conditions. One has to keep in mind that the realistic boundary condition depends on a variety of factors: the 
orientation of the edge, deformation of the chemical bonds near the edge, edge reconstruction, and possible chemical 
functionalization of the unsaturated bonds. Theoretical studies of these conditions were performed in several papers 

The physics of electrons near the armchair edge is simple: the edge always acts as a reflector of the incident 
electron current. The scattering is affected by details of the edge structure, such as C-C bond lengths near the edge 
and non-carbon radicals attached to the edge. Some additional details are provided in Appendix IB II 

The physics of Klein and zigzag edges, however, is quite different. These edges bind electrons. When the nearest- 
neighbor hopping Hamiltonian is used to describe graphene, the bound eigenstates (edge states) form a dispersionless 
band at the zero of energy (se e Appendix IB 2[) . These edge states can be observed experimentally as a peak in the 
local density of states [4J, [EH, [54[. For example, Fig. [5] shows scanning tunneling microscopy data from Ref. [54^ . 
There, the edge states are seen near the zigzag edge as a stripe of bright spots extending along the edge. 

The dispersionless band of the bound states is unstable with respect to different perturbations of H. For example, 
the inclusion of longer-range hopping makes the band disperse (Refs. foil l52l lo5l - l58| and Appendix IB 21) . 
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(a) 




Vs (V) 

FIG. 5: (Color online) Scanning tunneling microscope image of different graphene terminations, from Ref. (a) The 

termination of a graphene sample is investigated using a scanning tunneling microscope. Fragments of both zigzag and armchair 
types are identified. The edge states near the zigzag edge are clearly visible as stripes of bright spots stretching along the zigzag 
edge. No edge states are present near the armchair termination, (b) A typical dependence of the differential conductance near 
the zigzag edge is plotted. The peak near the zero voltage (V s = 0) corresponds to the edge states. Reprinted figure with 
permission from Y. Kobayashi, K.-i. Fukui, T. Enoki, K. Kusakabe, and Y. Kaburagi, Phys. Rev. B 71, 193406 (2005). 
Copyright (2005) by the American Physical Society. 

The most interesting way of lifting the degeneracy of edge states is by adding electron-electron interactions to 
H . It was predicted quite some time ago [59[ that magnetic correlations develop at a zigzag edge as a result of the 
interaction. This effect was investigated in several papers |60h63|| 

For example, a detailed DFT study was reported in [60|. It predicts that an isolated graphene zigzag edge is a 
ferromagnet with magnetic moment m of 0.3 of the Bohr magneton per unit cell of the zigzag edge. This may be 
understood qualitatively as follows. The spin coupling between nearest-neighbor carbon atoms is antifcrromagnctic. 
Thus, different sublattices have opposite magnetic momenta. In the bulk, this would lead to cancellation of the total 
moment. Near the edge, however, the electron density at the exposed row is higher that at rows located deeper 
into the bulk. Most importantly, all sites at the zigzag edge belong to the same sublattice; therefore, they have 
the same magnetic moment. This leads to a local imbalance of the total magnetic moment, which is seen as edge 
fcrromagnetism. 

Reference [60T ] reports many properties of the edge ferromagnetism (see also Table IIV[) : spin- wave dispersion 
(E = nq 2 , where k = 320meVA 2 ), stiffness (D = 2n/m = 2100 me V A 2 ), magnetic anisotropy (~ 10~ 4 ), and 
crossover temperature between the Heisenberg and the Ising regimes (T x ~ 10 K). The spin correlation length at room 
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Zigzag edge magnetism 


M, 

per unit cell 


T = 300 K 


T < 10 K 


Anisotropy 


0.3 /iB 


~ 1 nm 


~ 1 |iim 


Ising, 10~ 4 



TABLE IV: Summary of the magnetic properties of the graphene zigzag edge, as reported in Ref. [6Cj] . When the electron- 
electron interaction is taken into consideration, the edge-state degeneracy is lifted through the magnetization of the electrons 
near the edge. This creates a one-dimensional magnetic system. The magnetic momentum M of such system is 0.3 of the Bohr's 
magneton per zigzag edge unit cell. The correlation length f at room temperature is rather short, suggesting that it would 
be difficult to utilize the pristine zigzag edge in a spintronic device operating at room temperature. However, below T x — 10 
K, a crossover to Ising-like magnetic correlations occurs, and the correlation length increases exponentially upon approaching 
T — 0. It was proposed [13] that £ could be as large as a micrometer. 




FIG. 6: Graphene/graphane interface, as studied in Ref. [6(|. Larger dark balls correspond to carbon atoms, while the smaller 
light balls correspond to hydrogen. The lower part of the sample that is shown is graphene, while the higher part is graphane. 
In bulk graphane every carbon atom has a hydrogen atom attached to it. 



temperature is estimated to be of the order of one nanometer. At temperatures below T x it increases exponentially 
as the temperature decreases. The effects of edge disorder on the magnetic properties of the zigzag edge were also 
investigated. 

In Ref. (6l| magnetic properties of small fragments of zigzag edge were studied using the Hubbard model. Such 
model is relevant for systems with rough edges, consisting of alternating fragments of different terminations. It is 
demonstrated that a very short, of the order of three lattice constants, zigzag sequence is sufficient to generate a local 
magnetic moment. 

Some other interacting effects have also been studied. For exam ple, the influence of the long-range Coulomb 
interaction and doping on the edge magnetism is discussed in Ref. [63( using the Hubbard model with Coulomb 
interactions. The interaction of edge states with phonons was studied in Ref. [641 ]. 



C. Graphene/graphane interface 

Graphane [13, Hf| is a hydrogenated sheet of graphene. Unlike graphene, graphane is a semiconductor with a gap 
of the order of few eVs. The graphene/graphane interface, shown in Fig. [BJ can be viewed as a type of graphene 
edge: low-lying electron states in graphene decay exponentially inside the gapped media of graphane. As demon- 
strated by molecular dynamics simulations [r36| . the interface remains almost atomically sharp even at sufficiently 
high temperatures, which is an extremely attractive feature since it reduces scattering and simplifies the theoretical 
description. 

Depending on the orientation of the interface relative to the crystallographic axis of graphene, one can distinguish 
a zigzag- type interface (as in Fig. [5]), or an armchair-type (at the right angle to the zigzag). The zigzag interface 



supports edge states whose electronic and magnetic properties were investigated in Ref. [67 
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Besides hydrogenation, graphene may be subjected to fluorination in order to produce fluoridated graphene [68l472| . 
Like graphene, the latter is a semiconductor with a gap of the order of a few eV. A similar conversion occurs upon 
functionalization of graphene by nitrophenyl j73j . The properties of the interface between pure graphene and the 
functionalized material must be similar to the properties of the graphene/graphane interface. 

D. Fabrication of high-quality edges 

Most of the theoretical work so far has assumed that the edges of the nanostructures are atomically perfect. Needless 
to say, this is not easy to realize experimentally. However, recently, substantial progress in the area of high-quality 
edge fabrication has been achieved (e.g., [74l. l75j). 

In Ref. [zl[ a chemical method of deriving narrow graphene stripes with sharp edges was reported. A graphene 
sample was placed in a solvent and subjected to sonification. Strips with sharp edges and widths varying from 50 nm 
to sub-10 nm were extracted from the solution. The strips produced were used to fabricate a field-effect transistor-like 
device. 

In Ref. [lH it was experimentally demonstrated that during Joule heating of the graphene sample with disordered 
edges carbon atoms at the edge were vaporized, and sharp edges were stabilized. Model calculations shown that 
the edge defects were healed through point defect annealing and edge reconstruction. This process was modelled in 
Ref. |7g]. These findings suggest that many theoretical predictions dependent on the edge quality could be tested 
experimentally. 

IV. GRAPHENE NANORIBBONS 

Nanoribbons, which are strips of graphene, are among the most studied mesoscopic graphene structures. There are 
several reasons for this. First, they demonstrate unusual physical properties, for example, edge states, which might be 
used in future spintronics applications. Second, nanoribbons are easy to produce and demonstrate excellent transport 
properties. 

Third, they have an energy gap in their single-electron spectrum. This gap is a consequence of the electron 
confinement, and it is inversely proportional to the width of a nanoribbon [19l |20j. This suggests that, at least in 
principle, a nanoribbon with a desired value of the gap may be fabricated. Finite gap and high mobility are both 
very useful for the design of field-effect transistors (FET): they allow for large on/off ratios, small losses, and high 
operating frequencies. For example, a nanoribbon-based FET realized in Ref. [77| demonstrated an on/off ratio of 
about 10 6 at room temperature. Other nanoribbon- FET devices were described in Ref. [78|, [79|]. (Carbon nanotubes 
also have a gapped spectrum. However, their fabrication process is much more involved). 

Nanoribbons are usually classified by their type of edge; for instance, there are zigzag and armchair nanoribbons. 
Nanoribbons may also have disordered Q or more complicated regular types of edges (8G |. 

A. Zigzag nanoribbons 

Many interesting properties of zigzag nanoribbons are related to the presence of edge states in the nanoribbon 
electron spectrum. These states may be derived, together with other low-lying states, with the help of Eq. ^ 
supplemented by a boundary condition suitable for the zigzag edge [8l| . The resultant spectrum is shown in Fig. [7^,. 
Two almost dispersionless branches connecting the Dirac cones, K and K' , correspond to the edge states. They are 
analogous to the edge states discussed in Sec. IIIII 

Similar conclusions about edge states may be reached usin g first-p rinciple calculations. These reveal that the zigzag 
nanoribbon is a semiconductor with a width-dependent gap |82h84| . The lowest energy states are edge states. 

1. Edge magnetism 

Edge states are responsible for magnetism in zigzag nanoribbons. Edge magnetism is an interesting feature with 
potential spintronic applications; since the edge-state branch is both magnetized and able to carry current, it can be 
used to couple spin magnetization and current. This property may be used to control the magnetization with current 
or vice versa. 

The magnetism of a nanoribbon with pristine zigzag edges was investigated theoretically in 1996 j9lj . It is quite 
similar to the magnetism of an isolated zigzag edge: each edge of the nanoribbon has a finite magnetization, which is 



10 




-1.0 -0.5 0.0 0.5 1.0 

k (7i/a ) 

FIG. 7: (Color online) Single-electron energy spectrum of a nanoribbon calculated with the help of the Weyl-Dirac equation ([3|, 
from Ref. [8l| . Due to the finite width of a nanoribbon, the transverse momentum is quantized. Therefore, the nanoribbon's 
spectrum consists of a number of branches corresponding to different values of the quantized momentum. Since the graphene 
lattice is anisotropic, the spectrum of a zigzag nanoribbon [panel (a)] differs in several respects from the spectrum of an armchair 
nanoribbon [panels (b) and (c)]. For the zigzag nanoribbon the remnants of two Dirac cones, K and K' , are visible. The almost 
flat branch connecting K and K' corresponds to the edge states. It acquires a very weak dispersion due to interference of the 
edge states (exponentially) localized at the opposite edges of the nanoribbon. For the armchair nanoribbon both K and K' 
'coalesce' together. According to the calculations of Ref. [8l[, the armchair nanoribbon may be either a semiconductor with 
a small gap [panel (b)], or a gapless metal [panel (c)]. The size of the gap depends on the nanoribbon's width. However, 
more elaborate treatments accounting for the electron-electron interaction [85l. l8a|. electron-lattice interaction [iH. l87l |88|| . or 
longer-range hopping [89L |90| | proved that the gap is always non-zero. Reprinted figure with permission from L. Brey and 
H. Fertig, Phys. Rev. B 73, 235411 (2006). Copyright (2006) by the American Physical Society. 



induced due to the instability of a nearly-flat edge-state band. There is a non-zero coupling between the magnetizations 



of the two edges. In Ref. [91| it was shown that this coupling is antiferromagnetic, i.e., the magnetization vectors 
at the opposite edges are antiparallel. Such result is easy to understand. Consider the following two statements: 
(i) repulsive interactions between electrons on a half-filled bipartite lattice induce an antiferromagnetic correlation 
between sublattices; (ii) in case of the zigzag nanoribbon, one of its edges always terminates in atoms of A sublattice, 
the opposite edge terminates in atoms of B. As a consequence of (i) and (ii), the local antiferromagnetic tendency is 
translated into weak inter-edge antiferromagnetic interactions. 

Since the pristine zigzag edge is likely to be chemically unstable [4l|, |42[ , it is therefore important to study nanorib- 
bons with non-carbon atoms or functional groups attached to the edges. The case of a zigzag nanoribbon with mono- 
hydrogenated edges was discussed in Ref. [83[ . Such nanoribbons support edge states and have an edge ferromagnetic 
moment with antiferromagnetic coupling between the edges. This is consistent with the results of Ref. [4lj . where an 
extensive list of various edge types and their properties was presented. However, not all versions of functionalized or 
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FIG. 8: Zigzag nanoribbon with disparity between edges, from Ref. [§2]. The nanoribbon edges are parallel to y-axis. The 
left edge is monohydrogenated, while the right edge is dihydrogenated. Here, filled circles are carbon atoms, empty circles 
are hydrogen atoms. Since the opposite edges of a pristine zigzag nanoribbon have opposite magnetic momenta, a disparity 
between the two edges may induce a non-zero magnetization of the nanoribbon. Indeed, a local-spin-density approximation 
study [52] revealed that the nanoribbon in the figure possesses a finite magnetic moment. This result can be easily understood 
qualitatively. The monohydrogenated edge is ferromagnetic, while the dihydrogenated is non-magnetic. Thus, the whole system 
is ferromagnetic. Reprinted figure with permission from K. Kusakabe and M. Maruyama, Phys. Rev. B 67, 092406 (2003). 
Copyright (2003) by the American Physical Society. 



reconstructed zigzag edge support magnetism. 

A chemical way to produce a nanoribbon with finite magnetic moment was proposed in Ref. (92| : since the op- 
posite edges of the zigzag nanoribbon have opposite magnetic momenta, a disparity (e.g., a non-equivalent chemical 
functionalization) between the two edges may create a nanoribbon with non-zero magnetization. Local-spin-density 
calculations reported in Ref. [§2| proved that fact for the nanoribbon whose one edge is monohydrogenated, while the 
other is the dihydrogenated (see Fig. [5]). This result can be easil y u nderstood qualitatively. The monohydrogenated 
edge is ferro mag netic, while the dihydrogenated is non-magnetic [411 ]. Thus, the whole system is ferromagnetic. 

Reference [93| presented a very detailed DFT study of the effect the edge functionalization exerts on the zigzag 
nanoribbon's magnetism. The main focus there was the monohydrogenated zigzag nanoribbon, where some of the 
hydrogen atoms were replaced by other radicals. When the edges have non- identical chemical structure (i.e., one edge 
is purely hydrogenated, the other has some of its hydrogens substituted) , it was determined that a finite magnetization 
may be generated. The effect is particularly strong for the oxygen substitution. This result is consistent with Ref. [92j ■ 
It was also reported that such nanoribbon has different band gaps for different spin orientations. This state can be 
described as a spin-selective semiconductor. 

In addition, it was shown in Ref. [§3| that not only the magnetic properties but the band gap of a zigzag nanoribbon 
is sensitive to the chemical functionalization of the edges as well. For example, doping with oxygen may close the 
gap, provided that its concentration is sufficiently high. A variety of other effects dependent on the edge chemistry 
was also discussed. 



2. Half-metallicity 



A half-metal is a conductor whose charge carriers are fully spin-polarized. This is a very desirable property with 
potential applications to spintronics because it can be used to create a fully spin-polarized current. Reference (82j 
suggested that the application of a transverse electric field to a zigzag nanoribbon closes the gap for one spin orien- 
tation. The gap for another spin orientation is increased even more (see Fig. [5]). A similar conclusion was reached in 



Ref. |94j. Note that this idea is also based on inducing a disparity between the two edges: this time the transverse 



field is the agent producing the disparity. 
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FIG. 9: (Color online) Zigzag nanoribbon without (a) and with (b) a transverse electric field, from Ref. [82|. Symbols a (/3), 
in red (blue), represent the spin-up (spin-down) orientations of the edge band electrons. In the lower part of panel (a) both 
orientations are present: spin-up is on the left edge (L) of the nanoribbon, and spin-down is on the right edge (R). The letter 
'M' stands for 'middle' of the nanoribbon. In the upper part of (a) the energies of different edge states are plotted. On the left 
edge, the spin- up states are filled (E < 0), and the spin-down are empty. On the right edge, the situation is reversed. Panel 
(b) shows what happens when a transverse electric field is applied. The electrostatic potential pushes down the states on the 
left edge, and pushes up the states on the right edge. As a result, the density of states at the Fermi energy is zero (finite) 
for spin-up (spin-down) electrons. The system becomes half-metal: it is a conductor (metal), since there is a finite density of 
current-carrying states at the Fermi energy, yet, these states correspond only to one spin polarization (spin-down). Reprinted 
by permission from Macmillan Publishers Ltd: INaturel Y.-W. Son, M. L. Cohen, and S. G. Louie, Nature 444, 347 (2006), 
copyright 2006. 



The proposal [82j summarized in the previous paragraph was disputed in Ref. [95| . where it was argued that, when 
the transverse field is applied, two spin polarizations have different gap values. However, both gaps are finite for any 
value of the electric field, producing a spin-selective semiconductor, and not half-metal. Reference [95[ attributed the 
discrepancy to the artifacts of the computational technique of Ref. [82| . 

A suitable functionalization can enhance the half-metallic features of the zigzag nanoribbon [96j |. Also, Ref. [97j 
explored a wider range of chemical modifications of zigzag nanoribbons in search for robust half-metallicity. The 
results in Ref. [H, [97j suggested that the chemical modifications may by a powerful tool for the control of the zigzag 
nanoribbon's half-metallic properties. 



B. Armchair nanoribbons 



The electron properties of the armchair nanoribbon are simpler than those of the zigzag nanoribbon. Both tight- 
binding model and Weyl-Dirac equation calculations show that graphene armchair nanoribbons with pristine edges 
may be either in a semiconducting (finite gap, Fig. [7b), or metallic (zero gap, Fig. [7b), state with the gap oscillating as 
a function of the nanoribbon's width. A more accurate numerical study [18[ from 1997, which allowed for deformation 
of the carbon-carbon bonds dangling at the edges, proved that the metallic state is unstable: the dangling bonds 
deform, inducing a finite gap in the electron spectrum. Thus, according to Ref. fl8| . the armchair nanoribbon is 
always a semiconductor (see Fig. ITU1) . 

The above line of reasoning was generalized in Ref. [88j , where it was shown analytically that the metallic state of 
an armchair nanoribbon is generically unstable: the edge bond instability is only one possibility. It was also discussed 
[111 how the electron gap of a finite-length armchair nanoribbon can be effectively closed with the help of chemical 
modifications of the nanoribbon's edge. 

A gap may also be generated by electron-electron interactions [Hj], or longer-range hoppi ng [89l. r90| . In several 
papers, the electronic gap was determined with the help of first-principles techniques. Refs. [80, [87| reported DFT 
calculations of the gap for the armchair nanoribbon with different widths, with both pristine and monohydrogenated 
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FIG. 10: Armchair nanoribbon's gap as a function of the nanoribbon's width, from Ref. [l8l ]. The tight-binding calculation 
(dashed line) predicts that the gap of the armchair nanoribbon vanishes periodically as a function of the nanoribbon's width. 
However, the numerical calculations within a more elaborate model, which allows for deformation of the carbon-carbon bonds, 
prove that the gap is non-zero (albeit very close to it) for any width (solid line), although, the dependence on the width remains 
oscillatory. An analytical demonstration of the instability of the zero-gap state is given in Ref. (8£|. This instability occurs 
because the increase of the elastic energy due to the lattice deformation is smaller than the decrease of the electron kinetic 
energy due to the gap opening. Figure is reprinted from: M. Fujita, M. Igami, and K. Nakada, J. Phys. Soc. Jpn. 66, 1864 
(1997). 



Armchair nanoribbon energy gap 


Energy gap [eV] 


Bulk semiconductors with similar values of the gap 


Width [nm] 


Width [V3a ] 


0.7 


Ge, InN 


2-3 


8-12 


from 1.1 to 1.4 


Si, InP, GaAs 


1-2 


4-8 



TABLE V: Gap values for the armchair nanoribbons of different widths found using the density functional theory, from Ref. [80l ]. 
The quantity V3ao represents the lattice constant. The gap of the armchair nanoribbon oscillates with its width; however, 
the gap is never zero. The details of the oscillations depend on the edge functionalization (for example, in Ref. [80fl. both 
pristine and monohydrogenated armchair nanoribbons are discussed). The largest values of the gap are quite insensitive to 
functionalization. In order to have a gap value similar to known bulk semiconductors, a very narrow armchair nanoribbon must 
be used. These results were reconsidered in Ref. [13], where it was claimed that density functional theory underestimates the 
gap, and that the use of the so-called GW approximation [9S| is more appropriate. Within the GW framework, the values of 
the gap substantially increase. 



armchair edges. The results of Ref. [80] for the gap are summarized in Table [V] References [ljl H(| reported the 
experimental measurement of the ga p. T he experimental value for the gap was found to be consistent with the results 
of DFT calculations. However, Ref. (84| claimed that DFT underestimates the gap, and the use of the so-called GW 
approximation [98| is more appropriate. GW values of the gap are significantly higher than DFT values. 

The effect of the edge functionalization on the spectral gap was studied in Ref. (93|. The gap was found to be 
robust against functionalization. This is different from the case of zigzag nanoribbons whose gap is very sensitive to 
the chemical structure of the edges (see subsection IIV Al) . 

Thus, it is expected on the basis of theoretical studies that an armchair nanoribbon is a semiconductor with a 
width-dependent gap, whose value is rather insensitive to the edge chemical structure. 



C. Nanoroads 



So far we have assumed that graphene nanoribbons are formed by cutting a piece of graphene into a narrow strip. 
Another way to define a nanoribbon was proposed in Ref. [99| : to sculpture a graphene nanoribbon by removing 
hydrogen atoms along a narrow strip inside a wider graphane sample, as shown in Fig. llll In such a case, a nanoribbon, 
called nanoroad in Ref. j9{|, is bound by two graphene/graphane interfaces, which are discussed in subsection IIII CI 

As it was established in Ref. (66| . the graphene/graphane interface remains almost atomically sharp even at high 
temperatures. This makes nanoroads a promising candidate to observe ballistic transport. 

The magnetic and elect ronic properties of graphene nanoroads, including the effects of spin-orbit coupling, were 
discussed in Rcfs. In Ref. [99[ the electronic structure of nanoroads was studied using DFT. It was found 

that the armchair-type nanoroads are semiconducting. As for zigzag nanoribbons, when they are wide enough they 
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FIG. 11: (Color online) Nanoroads (i.e., graphene nanoribbons bounded by two graphene/graphane interfaces), from Ref. [9!|. 
Small dark balls are carbon atoms, while the large blue balls are hydrogen atoms. The armchair nanoroad is on panel (a), while 
the zigzag nanoroad is on panel (b). The quantities N a = 1, . . . , N and N z = 1, . . . , N characterize the width of the nanoroads. 
For the nanoroad on panel (a) N = 13, while N = 6 for the nanoroad on panel (b). Reprinted with permission from: K. Singh 
and B. I. Yakobson, Nano Lett. 9, 1540 (2009). Copyright 2009 American Chemical Society. 



demonstrate edge magnetism. Zigzag nanoribbons are semiconducting in the antifcrromagnctic state and metallic in 
the ferromagnetic state. 

Several other works also used DFT to analyze nanoroad properties. In Ref. [lOOf . the zigzag nanoroad stability and 
electronic structure were investigated, and it was established that even extremely narrow zigzag nanoroads are stable. 
Narrow nanoroads are always semiconducting due to P eierls instability, which opens a gap. A similar mechanism is 
responsible for the gap in polyacetylene. Reference [lOlj studied the adsorption of hydrogen on a graphene nanoribbon, 
also formulating the rules governing such adsorption. It was proposed to use such process in order to create narrow 
nanoroads. 

Reference [l02| showed that, due to enhanced spin-orbit coupling at the interface, a nanoroad might be used to 
convert spin polarization into valley polarization and vice versa. Such a device can operate at temperatures of about 
IK. 

The theoretical research summarized in this subsection indicates that nanoroads may be an attractive alternative 
to usual nanoribbons, able to sustain ballistic propagation of electrons, and also exhibit unusual spin features. 



D. Transport properties of nanoribbons 

For applications, such as FET [77l [78|. the transport properties of graphene nanoribbons must be investigated. A 
study of the conductance through a pris tine nanoribbon within the framework of the Landauer formalism was presented 
in Ref. jl03l | (see also the review [104| ). In such a case, the nanoribbon's conductance is quantized: it changes in 
discrete steps when the gate voltage is varied. For zigzag nanoribbons it was found that, when the gate potential is 
tuned to the charge neutrality point (i.e., the Fermi energy is at the apex of the Dirac cone), the conduct ance is finite 
due to the edge states. These are the only current-carrying modes under such conditions. Reference |103| | studied 
the transport through metallic (zero-gap) armchair nanoribbons. Since it is understood now that, strictly speaking, 
all armchair nanoribbons are semiconducting (see subsection IIV B| ) . the results obtained for such objects are valid as 
long as the gap may be neglected, e.g., when the temperature exceeds th e gap . Another study of electron transport 
through a disorder- free short -and-wide nanoribbon was pres ented in Ref. 1 1 051 ] . Its findings were compared well with 
the experiments in Ref. |l06| . It was concluded in Ref. |106| that the electron propagation in graphene is ballistic up 
to lengths of the order of 1 /im. 
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However, in a t ypical exp erimental situation ballistic propagation can be spoiled both by edge disorder |l07l4lld | 
and bulk disorder |l07lllld |. For bulk disorder, it was found that electron transport is rather insensitive to long-range 
disorder. Yet, when the disor der becomes short-range, it leads to Anderson localization and to destruction of the 
ballistic propagation |l07l . lllC| . 

If a nanoribbon is sufficiently narrow, the edge disord er is km Dortant [T07l - flTcill The armchair nanoribbons are more 
sensitive to edge disorder than the zigzag nanoribbons jl07l llldj . Reference [l08j | reported that when edge disorder is 
present, the difference between transport properties of the armchair and zigzag nanoribbons disappears. These results 
suggest that to e xperimentally produce a ballistic nanoribbon one has to overcome very stringent limitations on the 
edge purity [llCf . However, for a semiconducting armcha ir na noribbon of finite length one can optimize the width so 



that the l ocaliz ation effects are, to some exte nt, m asked [10 

Doping jlll| and edge functionalization (88l . lll2j affect the transport as well. For example, if a finite-length armchair 
nanoribbon has a gap induced by the edge-bond deformations (see subsection IIV Bp . then a suitably chosen chemical 
disorder at the edges may actually increase the conductivity by effectively closing the gap [H||. Namely, if the edges 
are functionalized with two different kinds of radicals randomly distributed along the length of the ribbon, then the 
term of the Hamiltonian responsible for the opening of the gap becomes disordered, and the gap closes when this 
term vanishes on average. Clearly, one has to counteract the Anderson localization in such a situation. Fortunately, 
this phenomenon is not important for a nanoribbon of sufficiently short length. Transport through this nanoribbon 
is effectively metallic |88| . 

Reference (86j pointed out that for nanoribbons with very corrugated edges the interaction effects can seriously 
affect the charge transport. According to Ref. [86|, a nanoribbon with corrugated edges can be viewed as a series of 
weakly coupled quantum dots defined by the random geometry of the edges. Electron transport through such system 
is limited by the Coulo mb block ade effect in these dots. 

Recent experiments 113l4ll6| on gated nanoribbons discovered that there is a Fermi energy interval wher e the 
conductance is suppressed, see Fig. [12j This phenomenon is called the transport gap. It was pointed out in Ref. |l!4{ 
that the size of the experimentally observed transport gap is too big to be consistent with conclusions of simple 
one-particle nanoribbon models, where the gap is generated due to the transverse quantization. To develop a more 
realistic description, the effect of both disorder and interactions must be accounted, in addition to the transverse 
quantization. It was proposed that the transport occurs through tunneling between consecutive "charge puddles" , 
which are the areas with non-zero charge induced by external the disorder potential. 

Our discussion shows that transport properties of the nanoribbons are affected most prominently by bulk and edge 
disorder, transverse quantization, edge type, and interactions. 



V. QUANTUM DOTS 

Quantum dots formed in semiconductor heterostructures have been studied ext ensively b ecause they are considered 
promising candidates for applications in optoelectronics on the nanometer scale [l 171 - 4121 . For instance, dots might 



be used in detectors, diodes, memory and laser devices. Furthermore, single-electron transport devices which make 
use of quantum dots could be employed as transistors, and spin-based dot devices might be useful for quantum logic 
gates. Electrons confined in usual semiconductor dots, with a typical size of a few hundreds of nanometers, are 
described by the Schrodinger equation and most of their electronic properties are now well-understood and have been 
experimentally studied by many research groups. 

The physics of graphene quantum dots is very different from that in usual semiconductor dots. The reason is 
twofold: (i) charge carriers in graphene are massless and obey the relativistic 2D Weyl-Dirac equation ([3]), and (ii) 
the different configurations of the carbon atoms at the boundaries of the dot affect significantly the dot properties. 
There are two basic methods of defining a graphene quantum dot. In the first method, the dots are defined by the 
actual geometry of the graphene layer and they are usually referred to as graphene islands. In the second method, the 
dots are defined through the application of electric and magnetic fields. Of course quantum dots can also be defined 
by combining these two methods and recently some other ideas have been put forward for dot formation, which for 
example include the application of strain to the graphene sheet, a spectral gap opening, and even chemical techniques. 



A. Geometry-induced dots and graphene islands 

It is now possible to mechanically cut (i.e., etch) a graphene flake into various shapes of a few tens of nanometers, 
which can confine electrons and thus act as quantum dots. These geometry-induced dots or graphene islands have 
well-defined discrete energy levels whose spectrum depends on the size, shape and the edge type of the dot. Further, 
disorder and interaction effects are also important for the electronic properties of any realistic graphene system. 
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back gate voltage (V) 

FIG. 12: (Color online) Graphene nanoribbon (a) and the nanoribbon's electrical conductance data (b,c), from Ref. |ll5l |. 
Panel (a) shows a nanoribbon defined by etching a graphene sheet. The darker area is graphene, while the lighter areas inside 
the red dashed lines are graphene- free. In panel (b) the differential conductivity is plotted as a function of the back-gate 
voltage and source-drain voltage. Note that for small source-drain voltage the conductivity remains suppressed (dark area) for 
small back-gate voltages, and grows when the latter exceeds a certain value. This is a manifestation of the transport gap. In 
panel (c) the same data is plotted in a smaller back-gate voltage window. The plot has the characteristic shape of overlapping 
noisy "Coulomb diamonds", suggesting that the transport occurs through several "charge puddles" acting as quantum dots. 
Reprinted with permission from: K. Todd, H. T. Chou, S. Amasha, and D. Goldhaber-Gordon, Nano Lett. 9, 416 (2009). 
Copyright 2009 American Chemical Society. 



A range of typical dot geometries including triangular, hexagonal , rectangular and circular have been studied nu- 
mericall y, mainly within the tight-binding and DFT models |l22l4l29| . In some cases exact analytical solutions are also 
possible [130l4l32| . For example, for a triangular armchair dot, exact tight-bi ndin g eigenfunctions and eigenenergies 
were obtained, and a technique for matrix elem ent calculation was developed [132| (see Fig. II 3|) . Quantum dots with 
arbitrary shapes have also been examined |l33| . 

An important feature of the nano-islands is the appearance of degenerate zero-energy states that are most l y loc alized 
at the edges, as predicted for triangular and circular zigzag dots, as well as rectangular dots 122Ml26l Il29l |. For 
triangular dots there is a sublattice imbalance, i.e., Nz=Na — Nsy^O, where Na (Nb) is t he nu mber of carbon atoms 
of sublattice A (B) , and this condition is sufficient in order to have Nz zero-energy states [l24| . The number of these 
states is proportional to the size of the edges which, in principle, can be made quite large. 

Nanostructures with degenerate zero-energy states are useful for applications, since the electrons inside such struc- 
tures may order magnetically. Magnetism is a consequence of the Coulomb interaction and Hund's rule. For example, 
the ground state of rectangular dots can support antiferromagnetic ordering whereby the magnetic moments are lo- 
calis ed at the zigzag edges with opposite orientation (for edge magnetism, see also the previous section). As shown in 
Ref. [l25j | for rectangular dots there is a critical minimum width between the zigzag edges that gives rise to magnetic 
ordering. If the width is smaller, then the state is nonmagnetic. On the other hand, triangular zigzag dots favor 
ferromagnetic ordering (see Table fVH . 

Interestingly, external uniaxial strain on square dots enhances the mag netization and leads to a spatial displacement 
(drift) of the magnetization from the zigzag to the armchair edges [123J. A magnetization enhanceme nt of 100% was 
predicted for a strain on the order of 20%, which might be possible to induce by mechanical methods |l34| . 

Moreover, it was theoretically predicted that the edge-state magnetism is robust to impurities a nd ed ge- defects, 
and survives even to irregular structures as long as there are three to four repeat units of zigzag edges (126| . A similar 
conclusion was reached in Ref. |6lj . Some examples of magnetic structures f or di fferent dots are presented in Fig. 1141 
The results are derived from the mean field analysis of the Hubbard model [l26j . 

Another promising property, which may be used in memory devices, is that the spin relaxation time can be long 
enough, as shown for triangular dots. In general, the spin relaxa tion time increases with the interaction strength and 
system size, but r emains long even for a relatively small system [l22| . 

Reference jl35| reported low-temperature electrical transport measurements on gated quasi-2D graphite quantum 
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FIG. 13: (Color onlin e) Th e Schrodinger equation for the tight-binding Hamiltonian on a triangular armchair quantum dot 
can be solved exactly [l32l ]. Moreover, the algebraic structure of the wave functions found is sufficiently simple to allow for 
analytical expressions for some matrix elements. In this figure the exact probability density for certain electron eigenstate on 
a triangular quantum dot is plotted. 



Graphene dots or islands (theoretical studies) 


Type 


Zero-energy edge states 


Magnetic ordering 


Triangular 


Yes 


Ferromagnetic 


Hexagonal 


No 


No 


Parallelogram 


No 


No 


Rectangular 


Yes 


Antiferromagnetic 



TABLE VI: Graphene dots have attracted considerable theoretical interest. Various geometries were examined and degenerate 
zero-energy edge states were predicted for triangular zigzag and rectangular dots (the rectangular dot is defined by two zigzag 
and two armchair edges). These states are mainly localized at the zigzag edges having only a small amplitude at the centre 
of the dot. Edge states are absent at arcmhair edges. The existence of such states leads to magnetic ordering that critically 
depends on the specific geo metr y. For rectangular dots there exists a minimum width between the two zigzag edges for stable 
antiferromagnetic ordering |l25l ] . The magnetic prop erties of the geometry-induced dots are robust to defects a nd i mpurities 
of the edges, can survive to irregular structures 126], and can be tuned by the application of an external strain |l23l |. 



dots. These were the first measurements on mesoscopic samples of graphite which consists of many stacked layers 
of graphene held together by weak van der Waals forces. Coulomb charging phenomena were demonstrated with 
the help of data in Fig. [T5l where the electrical current through the dot as a function of gate voltage and source- 
drain bias is plot ted. Mor e recent experiments probed the energy spectrum of quantum dots formed in a single 
layer of graphene 136l -ll38|. An all-graphene single- elect ron transistor, exhibiting Coulomb-blockade behaviour, was 



operational well-above liquid- helium temperatures [136j . The Coulomb-blockade peaks are (nearly) periodic as a 
function of gate voltage for large islands (> 100 nm), and nonperiodic for small ones (< 100 nm). The distance 
between the peaks is proportional to the sum of charging an d con finement energies. The former, being typically 
constant for a specific dot geometry, dominates for large islands jl36l | . For small islands the size quantization becomes 
important, and the confinement energy prevails, leading to nonperiodic peaks (see Table IVlT| . The energy-level 
statistics of g raphene islands was also probed, and it was shown to agree well with the theory of chaotic Dirac 



billiards [136 1 



Coulomb-blockade measurements on a graphene island (~ 200 nm) with an integrated charge detector were also 
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FIG. 14: (Color online) Magnetic properties of graphene dots were investigated numerically in Ref. |l26l ] within the framework 
of a mean field theory of the Hubbard model. Expectation value of the spin magnetization S z in four different dots. The size 
of the arrows is proportional to the magnitude of the magnetic moment. In this figure, N, N m and N„ denote the size of the 
dot in the perfect hexagon configuration, the number of armchair edges per vertical side, and the number of armchair edges 
per slanted edge, respectively. The magnetic moments are much larger at the zigzag edges than at the armchair edges and 
the internal sites. Reprinted with permission from S. Bhowmick and V. B. Shenoy, Journal of Chemical Physics 128, 244717 
(2008). Copyright 2008 American Institute of Physics. 



reported [139j. A nanoribbon placed 60 nm from the island acts as a detector, which enhanced the resolution of single 
charging events on the island. In addition, tunable double quantum dots w ere f abricated whereby the coupling to 
the leads and the interdot co uplin g were tuned by graphene in-plane gates |l40j . Spin spectroscopy has also been 
investigated in graphene dots jl4ll | . 



B. Field-induced dots 



Charge confinement within the "bulk" graphene sheet is tricky due to the Klein tunneling effect [142J. In case of 
normal incidence, this allows for perfect transmission of massless relativistic particles through high and wide potential 
barriers. A key point is that in the barrier region the states of massless Weyl-Dirac particles have an oscillatory 
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FIG. 15: (Color onlin e) Gated graphite quantum dots were fabricated and low-temperature electrical transport measurements 
were performed [l35l |. (a) Current versus gate voltage V g with source-drain bias V s a = 10 ^iV at temperature T ~ 100 mK. 
Coulomb oscillations are observed with a period in gate voltage of AV g = 1.5 mV. (b) The differential conductance (dI/dV s d) 
is plotted as a color scale versus gate voltage (V 9 ) and source-drain bias V a d- Blue (red) signifies low (high) conductance. The 
charging energy of the dot is equal to the maximum height of the diamonds: AV s d = 0.06 mV. Reprinted with permission 
from J. S. Bunch, Y. Yaish, M. Brink, K. Bolotin, and P. L. McEuen, Nano Letters 5, 287 (2005). Copyright 2005 American 
Chemical Society. 



Graphene dots or islands (experimental studies) 


Size of island (D) 


Coulomb peaks 


Energy scale 


< 100 nm 
> 100 nm 


Nonperiodic 
Periodic 


Confinement ~ VFh/2D w 41 meV (D = 50 nm) 
Charging « 3 meV (D ~ 250 nm) 



TABLE VII: Graphene islan ds we re investigated experimentally via electrical transport measurements, and single-electron 
transport was demonstrated 136]. When the diameter of the island is large (D > 100 nm) the Coulomb peaks in the 
conductance, as a function of back-gate voltage, are periodic and their position is determined mainly by the characteristic 
charging energy. For small-diameter islands (D < 100 nm), the position of the peaks is nonperiodic and the dominant 
energy scale is the confinement energy on the order of Ed ~ vph/2D (h is Planck's constant). This is much larger than the 
corresponding energy Es ~ h 2 /8mD 2 of Schrodinger electrons with effective mass m; Ed/Es ~ 40 for D = 100 nm, and m is 
the effective mass for GaAs. Stable, robust and conductive dot islands as small as 15 nm were fabricated, showing the potential 
of graphene for nanoelectronics. 



character, even at energies lower than the potential height, as happens exactly outside the barrier. This is completely 
different from Schrodinger particles with non-zero mass, for which the states in the barrier region decay exponentially 
and therefore perfect transmission is not feasible. Experimentally, the Klein tunneling was demons trate d in graphene 
through electrical transport measurements in steep potential barriers generated by metallic gates [l43j . 

Because of the Klein tunneling, an electrostatic potential minimum in graphene leads to quasi-bound states, i.e., 
reson ant states (see Table IVIIII for a classification of the dot states) and therefore it is inadequate to confine elec- 
trons jl44| - |l4d |. Nevertheless, the finite lifetime of the states, characterizing the trapping time of an electron in the 
dot region, can be relatively long. It depends on the potential profile and eigenenergy. A smooth potential and a 
large angular momentum enhance the lifetime of the quasi-bound states |144| . Il46| , as happens also with eigenenergies 
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Back gate V bg (V) Bias dot V b qd (mV) 

FIG. 16: (Color online) Graphene quantum dot device with charge detector and transport measurement data [1391 ] . (a) Scanning 
force micrograph of the measured device. The central island, that acts as the quantum dot (QD), is connected to source (S) 
and drain (D) contacts via two narrow constrictions. The diameter of the dot is 200 nm and the constrictions are 35 nm wide. 
The charge detector (CD) is a graphene nanoribbon and the lateral gates Bl, B2 and PG are used to tune the device, (b) 
Current as a function of back-gate voltage of the QD (upper panel) and CD (lower panel). The source-drain voltage is 500 /iV 
and the temperature is 1.7 K. The inset shows typical Coulomb-blockade features as expected in a dot device, (c) Differential 
conductance is plotted as a color scale versus source-drain bias and PG gate voltage, for a back-gate voltage of 2 V. The 
charging energy was estimated to be about 4.3 meV. Reprinted with permission from J. Guttinger, C. Stampfer, S. Hellmuller, 
F. Molitor, T. Ihn, and K. Ensslin, Applied Physics Letters 93, 212102 (2008). Copyright 2008 American Institute of Physics. 



close to the maximum of the potential barrier [l45| . Figure [17] shows the effect of the quasi-bound states on the local 
density of states of a circular dot. 

The electrostatic confinement of electrons in graphene dots was a lso examined through the dependence of the 
conductance on the dot area, which is tunable with a metal gate |l47| . Both disc-shaped dots, in which the classical 
dynamics is regular, and stadium-shaped dots where the classical dynamics is chaotic were studied. Confinement can 
be achieved only in the former when the corresponding Weyl-Dirac equation is separable. 

The application of a magnetic field can completely suppress the Klein tunneling, leading to bound states [l48l ll 50l — 

Il52j . Thus graphene quantum dots can be formed with the help of a nonuniform magne tic fi eld, whereby the field is 
zero within a disc area defining the spatial region of the dot, and nonzero outside the dot |l50j |. The combination of an 
electrostatic poten tial and a vector magnetic potential allows the confined-deconfincd character of the dot states to be 
tuned at will |l48l ] . Most interestingly, it allows gra phene dots to be formed in a uniform magnetic field using standard 
gate electrodes as in common semiconductors [1481 ] . Then the quantum states can be tuned with the strength of the 
magnetic field and this property allows the Klein tunneling me chan ism to be probed experimentally in graphene dots. 
A dot design suitable for this experiment wa s sug gested in Ref . |l48l | . The concept of defining a magnetic graphene dot 
was further developed theoretically in Ref. [l49(. In particular, in a strong magnetic field the electrostatic potential 
of the dot is adjusted so that the confined dot states lie in the gap between Landau level and Landau level —1 (see 
Fig. n~9|) . This ensures that the dot states are energetically isolated in a region of low density of states and thus they 
can be probed using standard charge-sensing measurements as in a GaAs dot. Numerical esti mate s showed that a 
typical spacing between the dot levels is ~ 2 meV at a magnetic field of 5 T. In addition, Ref. |l49j [ considered how 
this basic idea can be extended to a graphene antidot for which the levels of the confined states lie in the region 
between Landau level and Landau level +1. For a confined state with energy E in the dot, there is a corresponding 
confined state with energy — E in the antidot. This is a unique property of graphene due to the symmetry of the 



21 




FIG. 17: (Color onlin e) T he physics of quasi-bound states in circular graphene dots was examined theoretically by solving the 
Weyl-Dirac equation [l45l |. The figure shows the local density of states as a function of energy in the case of a barrier height 
of V = 12 (dimensionless units) for three angular momentum numbers: m = (solid curve), m — 1 (dotted curve), and m = 2 
(dashed curve). The peaks become narrower as the momentum increases, within a specific energy range, thus the lifetime of the 
corresponding states becomes longer. Notice the very narrow peak when the energy is close to the barrier height (for m = 0), 
as a consequence of the total reflection of the wavefunction at the dot edge. Reprinted figure with permission from A. Matulis 
and F. M. Peeters, Physical Review B 77, 115423 (2008). Copyright (2008) by the American Physical Society. 



Quantum dot states 


Type 
(other used names) 


Lifetime 


Spatial asymptotic behaviour of wavefunction 


Bound 
(confined, stable) 
Quasi- bound 
(deconfined, resonant) 


Infinite 
Finite 


Exponential decay 
Oscillatory 



TABLE VIII: The quantum states of a dot, formed within the "bulk" graphene sheet, can be either bound or quasi-bound. 
Because of the Klein tunneling, both types of states can have a large amplitude in the barrier region of the quantum dot, though 
their asymptotic behaviour is different. Exponential decay is characteristic of bound states, whereas oscillatory behaviour is 
characteristic of quasi-bound states. As shown in Ref. |l48l | the type of states can be tuned with an electrostatic potential and 
a uniform magnetic field perpendicular to the graphene sheet (see also Table llX| . 

Dirac cone. The physics of graphene antidots in a magnetic field was also examined in Ref. [l53j . 

Figure [T8l illustrates the magnetic field-induced suppression of the Klein tunneling for an electron inside a do t, and 
Table ITXl summarizes the general conditions for confinement in a circular quantum dot, as derived in Ref. |l48j |. The 
physics of electrostatic barriers in the presence of uniform and nonuniform magnetic fields is analyzed in Sec. I VIII 

C. Nanoribbons of graphene and dots 

As discussed in the previous section, the most commonly studied graphene nanoribbons are either of armchair 
or zigzag type. A quantum dot defined by an external electrostatic (parabolic) potential along the nanoribbon was 
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FIG. 18: (Color online) Klein tunneling in a circular graphene quantum dot in a uniform magnetic field. The dot is defined by 
the electrostatic potential V(r), which vanishes at the centre (O) of the dot and asymptotically rises to a constant value, while 
the magnetic field B is perpendicular to the graphene sheet, (a) For B = and small angular momentum, the radial probability 
distribution p(r), for one of the spinor components, has a large amplitude near the centre of the dot and oscillates inside the 
barrier region because of the Klein tunneling. In this case, which is unique to graphene, the quantum state is quasi-bound and 
has an oscillatory asymptotic character, (b) When the magnetic field is nonzero, the Klein tunneling is partially suppressed. 
At large r the oscillatory behavior is replaced by exponential decay, indicating that the state is bound. As in (a), this case is 
also unique to graphene. (c) With increasing magnetic field, the Klein tunneling is completely suppressed and the probability 
distribution decays exponentially inside the electrostatic barrier. The state is now bound near the centre of the dot and such a 
state can be seen in both graphene dots and usual semiconductor dots, e.g., GaAs. For the latter, this state can be seen even 
when B = 0. A magnetic-field- indu ced confinement-deconfinement transition in a graphene dot due to the Klein tunneling was 
theoretically examined in Ref. [l48| . 



investigated theoretically in Ref. (154| . Quasi-bound states inside nanoribbons with either zigzag or armchair edges, 
in both metallic and semiconducting samples, were predicted. The lifetime of such states can be made long enough 
by increasing the characteristic length of the external potential, for example, with a gate electrode. Further, the 
dependence of the conductance on the gate voltage was found to be sensitive to the type of edges as illustrated 
in Fig. [20] The proposed device can operate either as a quantum dot or as a point contact. For the former, the 
conductance displays resonances at negative Fermi energy, whereas for the latter it has a step-like behaviour at 
positive Fermi energy. The conductance steps in the semiconducting system are twice smaller than those in the 
metallic system. 

Se miconduct ing nanoribbons with armchair edges were proposed for the formation of spin qubits in graphene 
dots jl55l Il56| . Confinement in one direction is achieved naturally by the nanoribbon and in the second direction 
electrically by gate voltages. In this set-up, the valley degeneracy is lifted, thus allowing Heisenberg spin exchange 
coupling in tunnel-coupled dots. Such graphene dots can be coupled over long distances as a consequence of the 
rclativistic nature of electrons in graphene, exhibiting Klein tunneling. 

The electrostatic confineme nt o f electrons in graphene nanoribbons as well as the Coulomb-blockade effect were 
experimentally demonstrated [l57l |. In particular, electrons are confined between gate- induced pn-junctions acting as 
barriers. However, even when no pn-junctions are formed, the electrons are still confined, though in a larger area due 
to strong disorder. 
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FIG. 19: A grapheme dot or an antidot can be formed with the help of a uniform magnetic field applied perpendicular to the 
graphene sheet [149| | . (top) The electrostatic potential profile of an antidot (solid line) and a dot (dashed line) along the radial 
direction. Energy spectra as a function of the applied magnetic field of dot (left) and antidot (right). The dashed lines show 
the Landau levels of an ideal graphene sheet. The potential is adjusted so that the confined dot (antidot) states lie in the gap 
between Landau level and Landau level —1 (+1). The symmetry between the energy levels (E — > —E) of the two systems is 
a direct consequence of the Dirac cone band-structure in graphene. 



Circular graphene 


dot: V = V r% A e = A a r f 


Asymptotically 


s, t 


Dot states 


Ag < V 


t < s 


Quasi-bound for all Vo, Ao 


Ag~V 


t = s 


Quasi-bound for Vo > VFeAo 


Ag~V 


t = s 


Bound for Vo < VFeAo 


Ag>V 


t> s 


Bound for all Vo, Aq 



TABLE IX: Confinement of electrons in a circular graphene dot is conditional because of the Klein tunneling [l48t |. Consider a 
graphene dot defined by the electrostatic potential V = Vor s and the magnetic vector potential A = (0, Ag,0), with Ag = Aor 
and s, t > 0. In such a situation, if asymptotically Ag < V, then the quantum states are quasi-bound; if Ag > V, they are 
bound. In the special case Ag ~ V , the states are bound only when Vo < VFeAo (e is the absolute value of the electron charge). 
This suggests that the states can be tuned at will with an electric or magnetic field. A simple limit occurs for a constant 
potential (V = const.) and uniform magnetic field B perpendicular to the graphene sheet (t = 1, Ao = -j). Then the states 
are always bound and correspond to the well-known Landau levels. 
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FIG. 20: Theoretic al st udies showed that a graphene quantum dot can be formed by an external electrostatic potential in a 
nanoribbon system [l54l ]. The conductance through the dot as a function of the Fermi energy depends on whether the system is 
semiconducting (solid line) or metallic (dashed line). In both cases, the nanoribbons have armchair edges. For e < the device 
operates as a quantum dot and the conductance exhibits resonances. For e > the device operates as a point contact and the 
conductance exhibits steps. For e > the conductance steps in the semiconducting system are twice smaller than those in the 
metallic. Reprinted figure with permission from P. G. Silvestrov and K. B. Efetov, Physical Review Letters 98, 016802 (2007). 
Copyright (2007) by the American Physical Society. 



D. More dots 



Tunable quantum dots w hich take advantage of a gap in the energy dispersion were also proposed in both mono- 
layer and bilayer graphene |l58l4l61 |. Experimentally, the gap can be introduced via a chemical and/or electrical 
technique [ID, H|| (see also Table |TT| , and it allows dots to be formed electrostatically in a quite similar manner 
as in common semiconductors. A finite gap introduces a mass term in the Weyl-Dirac Hamiltonian. Then for a 
dot-confining potential with a finite asymptotic value, the gap gives rises to an energy range within which the Klein 
tunneling is suppressed, leading to the formation of bound states. In this energy range, which is directly proportional 
to the value of the gap, hole states do not exist and therefore the electron states decay exponentially. Moreover, for 
quantum dots formed in the gapped sample the valley degene racy is lifted by a uniform magnetic field. This property 
might be attractive in order to define spin and valley qubits 158 1. 

It was also shown theoretically that a spatially modulated Dirac gap in the graphene sheet can lead to confined states 
with discrete energy levels, thus giving rise to a dot. The basi c advantage of this proposal is that the dot is formed 
without applying external electric and/or magnetic fields jl6fj |. Thus magnetic fields can be used to manipulate the 
spin states without affecting the confinement of the corresponding orbital states. The properti es of a Dirac gap-induced 
graphene dot in the presence of an electrostatic quantum well potential were studied in Ref. |l6l| . It was shown that 
confined states which are induced thanks to the spatially modulated Dirac gap couple to the states induced by the 
potential. The resulting hybridised states are localised in a region which can be tuned with the potential strength; 
an effect which i nvolv es Klein tunneling. Numerical calculations of the local density of states suggest that this effect 
could be probed [l6l| . 

Strain engineering is another proposal in order to generate confinement in a sheet of graphene [162| . Local patterning 
of the substrate induces in-plane strain in the graphene lattice, anisotropically changing the hopping amplitude 
between the carbon atoms. As a result, in the continuum approximation the quasi-particles are described by an 
effective Weyl-Dirac equation in the presence of a gauge field. It turns out that this field can act in a rather similar 
manner a s a m agnetic field and therefore it can lead to confined quantum states. A noteworthy advantage of this 
proposal [162] is that patterning can be made directly on the substrate, hence protecting the graphene layer from 
possible damage. 

Vacancy clusters in the graphane sheet were also suggested for dot formation. In particular, DFT and tight binding 
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FIG. 21: (Color online) The Aharonov-Bohm (AB) effect was experimentally demonstrated in a graphene ring device [164 
(a)-(c) AB conductance oscillations versus applied magnetic field, for a back-gate voltage of 30 V and temperature of 150 mK. 
For B ~ 3 T, an increase of the AB amplitude is observed. Panel (a) shows a magnified view of the yellow region in (c), while 
panel (b) expands the blue part in (c). Reprinted figure with permission from S. Russo, J B Oostinga, D. Wehenkel, H. B. 
Heersche, S. S. Sobhani, L. M. K. Vandersypen, and A. F. Morpurgo, Physical Review B 77, 085413 (2008). Copyright (2008) 
by the American Physical Society. 



calculations showed that cluster of hydrogen vacancies can serve as quantum dots. The stability as well as the shap e 
and size of these dots depend crucially on the graphcne/graphanc interface energy and the degree of aromaticity [163| . 



E. Quantum rings 



Quantum rings in graphene have also attracted some interest, mainly because these types of devices allow the 
investigation of phase-coherence phenomena, as it is now well-known from studies in usual semiconductor systems. 

In a graphene ring d evice , conductance oscillations versus magnetic field were reported as a consequence of the 
Aharonov-Bohm effect |164| . The amplitude of the oscillations increases at high magnetic field in the regime where 
the cyclotron diameter becomes comparable to the width of the arms of the ring. For temperatures below 1 K the 
extracted phase-coherence length is comparable to or larger than the diameter of the ring, which is approximately 1 
/im. 

Theoretical investigations of graphene rings showed that the valley-induced orbital degeneracy is lifted, as a result 
of the ring confinement and the applied magnetic field 165]. This lifting has observable consequences on the persistent 
current and the ring conductance. An interesting finding is that the degeneracy can be controlled with the induced 
Aharonov-Bohm flux, and this can be achieved irrespective of the magnitude (weak or strong) of the intervalley 
scattering. 

Another theoretical wor k sh owed that both electrons and holes can be confined in electrostatically formed quantum 
rings in bilayer graphene |l66j . There are two main advantages in this proposal. First, bound states can be created 
owing to a position-dependent energy gap that suppresses the Klein tunneling. Second, the ring parameters can be 
tuned by external fields. 
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FIG. 22: (Color online) The pn-junction studied in Ref. [lT3t ] . There are two gate electrodes in this device. The first one is 
the semi-infinite gate on the left side. It controls the density drop p2 — pi across the junction. The second electrode is an 
infinite back gate above the sheet (not shown). It fixes the density p2 at far right. Lines with the arrows show trajectories 
of an electron (— ) and a hole (+). The electron current in n-region is converted into hole current in p-region. Note that the 
direction of the incident electron current in n-region and the direction of the hole current in p-region are symmetric with respect 
to j/-axis reflection (the same type of refraction is shown in Fig. [23}. Reprinted figure with permission from L.M. Zhang and 
M.M. Fogler, Phys. Rev. Lett. 100, 116804 (2008). Copyright (2008) by the American Physical Society. 



The role of Coulomb-induced electron-ele ctron interactions and their interplay with the valley polarization in a 
graphene quantum ring were also examined |l67l |. In a few-electron ring, the interactions have a direct signature on 
the fractional nature of the Aharonov-Bohm oscillations in the persistent current and the absorption spectrum, and 
therefore they could be observed. 



VI. GRAPHENE PJV-JUNCTIONS AND PNP-STRUCTURES 

Several graphene-based field-effect devices have been realized in laboratories. Reference [l3| reported t he fa brication 
of a FET made of graphene which operates at a record-breaking frequency of 100 GHz. In Ref. jl68| | a room- 
temperature-operated switch demonstrating an on/ off r atio exceeding 10 6 was described. Also, the implementation 
of a digital integrated circuit was reported in Ref. [l69j | . The microcircuit consists of two tran sistors and performs 
the logical inversion operation. A graphene FET used as a biosensor was described in Ref. |l70j . 

All such devices are characterized by a spatially inhomogeneous Fermi level inside the graphene sample. There are 
several basic types of such systems: interfaces separating regions with different concentrations of the charge carriers 
(pp' -junctions, nn'-junctions), or regions with carriers of opposite signs (pn-junctions), or series of such interfaces (pnp- 
structure, pp'p-structure, etc.). Transport properties of these systems is an important subject of both theoretical and 
experimental investigations. These studies are reviewed below. 



A. pn-j unction 

If two planar electrostatic gates separated from the graphene sample by an insulating layer are charged in such a way 
that the chemical potential at x > is shifted above the electroneutrality level and at x < below the electroneutrality 
level, a pn-junction i s formed at x — 0, see Fig. [521 The simplest model of the electrostatically-defined pn-junction 
was studied in Refs. [ml Ell- 
in Ref. |l7l| the current transmission through a pn-junction was investigated. When the current approaches the 
pn-junction at the right angle, it passes through with no reflection. This is a manifestation of the Klein tunneling. 
Otherwise, the reflected current appears. The primary source of reflection in such a system is a classically-forbidden 
strip near the center of the junction, which can be crossed only by quantum tunnelling. The strip's width I de pend s 
on the incidence angle 9 [for normal incidence 1(6 = 0) = 0, hence, Klein tunneling]. In the model of Ref. |l7l| . 
parameters of the potential barrier under which the particle has to tunnel depend on the geometry of the pn-junction 
and the gates' potentials. 

The main finding of Ref. jl7l| was that the current transmission through such pn-junction is very sensitive to the 
angle 9: for normal incidence the transmission is perfect, but it exponentially quickly deteriorates when 9 grows. This 
allows one to create very collimated beams of current. Additionally, the selectivity to 9 can be used to detect the 
magnetic field: since the magnetic field bends the trajectory of a charged particle, then, in a properly designed device, 
a particle hits the interface at a magnetic-field-dependent angle. As a result, the transmission b ecom es sensitive to 
the field. Different devic es ut ilizing properties of the gr aphe ne pn-junction were proposed in Ref. |l7lj . 

The treatment of Ref. [l7l| was re-examined in Ref. (173| | . It was noted there that the non-linear charge screening 
affects the pn-junction characteristics. When the non-linear screening is accounted, significant deviations from the 
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FIG. 23: (Color online) pn-junction as an electron-focusing device, as described in |l79l |. If the gates' voltages are such that 
the concentration of electrons in the n-region and the concentration of holes in the p-region are the same, a particle hitting the 
pn-junction interface at the incidence angle C would be converted into a hole on the other side of the junction propagating at 
the refraction angle V — —6 C - Under such condi tions , a current emanating from a point source is focused into a small spot, 
"the image", on the other side of the p?vjunction [l79l |. 



findings of Ref. jl7l| . which neglects many-body effects, are discovered. The results in Ref. |l73l | indicate that the 
interaction significantly reduces the pn-junction resistance. 

A mo re ge neral study, including not only the electron-electron interaction, but also the disorder, was presented 



in Ref. [17J]. It was shown that, depending on the junction's parameters, it may be in either of three regimes: 



(i) ballistic, where the pn-junction resistance is dominated by the ballistic contribution, (ii) diffusive, where the 
resistance is dominated by the diffusiv e con tribution, and (in) the crossover re gime, wh en both ballistic and diffusive 
contributions are comparable. In Ref. |l74l | several experimental pn-junctions |l75l - ll78j were analyzed trying to find 
junctions in the ballistic regime. It was concluded that the considered experimental systems satisfy the conditions for 
ballistic propagation only marginally at best. It was suggested that higher mobility or a larger carrier concentration 
gradient near the junction is required to create a ballistic pn-junction. 

T he ba llistic pn-junction has attracted considerable attention due to its unusual electron-refraction properties. In 
Ref. |l79j it was observed that, under certain conditions, the electron beam passing through a graphene pn-junction 
experiences refraction in such a manner that the refraction angle equals to minus the angle of incidence, see Fig. 1231 
If this is the case, then current emanating from a point source on one side of the pn-junction is focused into "a point 
image" on the other side of the junction. This situation is similar to the refraction of light at the interface with a 
mctamaterial whose refra ction index is minus unity. The focusing properties of the pn-junction with circular geometry 
were investigated in Ref. [l8dj . 

However, this ability to focus the electrical current is easy to spoil. Ref. jl79| pointed out that at the level of the 
geometrical optics the focus is p erfect only if the density of holes in the p-electrode is the same as the electron's density 



in the n-electrode. In Ref. |174j it was shown that disorder destroys the focus as well. Finally, since the transmission 



of a pn-junction decays quickly as the incidence angle deviates from 7r/2, only a small fraction of the initial current 
is able to pass through the junction to form "the image" . 

When a graphene pn-junction is placed in a non-uniform magnetic field, it acquires new interesting features. This 
type of devices is discussed in section IVHI 



B. Doping graphene by contact with metals 



In addition to electrostatic doping, it is possible to change the charge density in graphene by making contact with a 
metal electrode. In such a case, depending on the electrode's material, the electrons either leave the graphene sample 
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to the electrode or flow into the graphene from the electrode. 

Junctions created with the help of this kind of doping were investigated experimentally in Ref. [18 lj ] . Materials for 
the metallic electrodes were chosen in such a way as to dope the graphene with holes. Then, depending on the voltage 
of the back gate, either pp'- or p?i-junctions were formed. 

The graphene- metal interface was investigated theoret ically in Refs. |l82l4l84 |. Charge transfer between t he m etal 
electrode and the graphene sample was studied in Refs. |l85l . Il86l | with the help of DFT. According to Ref. |l85j | the 
Fermi energy shift inside the graphene sample is a monotonous function of the metal work-function, as one should 
expect. However, when the work-function of the metal coincides with that of graphene, the graphene sample is not 
neutral, as one naively might expect, but rather it is predicted that the sample is n-doped. This happens because of 
the chemical interaction between the metal and graphene. 



C. prep-structure 



The theory of electronic transport in clean pnp-structures was presented in Ref. [172J. There, they demonstrated 
that the conductance of ballistic jmp-structure exhibits oscillations ('Fabry-Perot' resonances) as a function of the 
carrier concentration in the middle (n) area of the pnp-structure. These resonances are due to quasi-bound electron 
states in the n-region of the pnp-structure. 

A more general numerical study, which accounts for interaction and impurities, was performed in Ref. jl87l |. It 
reported a crossover from ballistic to diffusive regime when the mean-free-path becomes comparable to the length of 
the middle region. The disorder wipes out the 'Fabry-Perot' resonances. However, it is conceivable that these survive 
under a small concentration of impurities, and, thus, could be seen experimentally. 

A phenomenon analogous to the 'Fabry-Perot' resonances was discussed in Ref. [l88l ll8P| . where the transmission 
through several junctions connected in series was studied. Because of electron wave function interference, the transport 
through such structure demonstrates a non-monotonous dependence on the current incidence angle and the distance 
between the junctions. 



Electrostatically-defined npn- and pnp-structures were realized experime ntall y |l90l - ll92' |. For example, Fig 



shows a scanning electron microscope image of a pnp-structure f rom Ref. [1901 ] . Th e exp erimental observation of 
'Fabry-Perot' oscillations in pnp-structures was rep orted in Ref. [ 1 9 OL 1 1 9 1| . In Ref. jl92| | experimental data were 
analyzed within the theoretical framework of Ref. j 1 T lj ] . Reference [l92t ] concluded that, in the fabricated pnp- 
structure, the individual pn-junctions are ballistic, and that the fabrication of a ballistic graphene pnp-structure is 
feasible. 

In the previous subsection we discussed the peculiar electron refraction at the pn-interface. In Ref. [l79| several 
possible applications of this effect were proposed, among which the most known is the so-called "electron Veselago 
lens" . The latter device is a ballistic graphene npn- or pnp-structure in which both junctions are tuned to operate in 
such a manner that the electronic current emitted from a point current source in the left n-electrode travels through 
two junctions and would be focu sed into a point (image of the source) in the right n-electrode. A similar phenomenon 
was predicted by Veselago [193| in the optics of materials with a negative-refractive index: the electromagnetic rays 
emitted from a point source are focused upon passing through a slab of such material. This slab is called the Veselago 
lens. It is an analog of the npn-structure under discussion. 

However, t he above analogy is inco mplete. The "superresolution" , the most a dvant ageous property of the optical 
Veselago lens 

HMDS! (see also Ref. fl96Hl99| '). is absent for the graphene device j20d| . Also, since the current trans- 
mission decays quickly for non-normal incidence, the graphene lens is very opaque. This might make its application 
problematic. 

Our discussion in this section suggests that graphene pn-junctions and pnp-structures, due to their interesting 
properties, may, in principle, be used for current control and magnetic field sensing applications, provided that a 
way to attenuate the effects of the disorder is found. If a graphene pnp-structure is placed into a non-homogeneous 
magnetic field, it may, under certain condition, act as a electron waveguide. The relevant discussion can be found in 
Sec. rvTTCl 



VII. QUANTUM BARRIERS, WIRES, AND WAVEGUIDES 



In this section, we address charge transport for designing tunable charge-conducting elements. Unlike quantum dots 
(Section V), where electrons are bound in a closed space, in graphene-bascd quantum wires and waveguides the charged 
particles should be confined only in one direction and be freely propagating in another one, as in pnp-structures. 

There are several methods of charge confinement in graphene (see Sections IV- VI ab ove) . Typically, this is achieved 
either (i) chemically (by binding graphene atoms to foreign atoms: e.g., oxygen j20lj . fluoride [68|, O, [72[ , hydrogen 



29 





FIG. 24: (Color online) Experimental realization of a graphene pnp-structure, from Ref. [1901 ] . Panel (a) shows the scanning 
electron microscope image of the structure. Large purple rectangle is the graphene sheet. Two bulk yellow electrodes (source 
and drain) and one narrow blue electrode (top gate) are placed on top of the graphene. The back gate beneath the structure 
is not visible. The inset presents an enlarged view of the top gate. Panel (b) shows a schematic diagram of the same setup. 
Reprinted by permission from Macmillan Publishers Ltd: Nature Physics, A. F. Young and P. Kim, Nat. Phys. 5, 222 (2009), 
copyright (2009). 



or aryl groups [73j); (ii) mechanically (either cutting or bending graphene sheets [162L 12031 12041 ] . or by 



creating inhomogencous spatial strain distributions [HI, |204| - |206( ); (Hi) thermally 12 1, and (iv) electronically (by 
applying electromagnetic fields). 

Each of these methods has its advantages and disadvantages. The first two are very effective; however, they 
are rather difficult to control, in the sense that any tuning (change of parameters) requires a reconstruction, either 
chemical or geometrical, of the whole graphene sample, which usually cannot be done quickly. Thus, these methods 
are not very suitable in designing tunable electronic devices. 

Here we concentrate on confining electrons using electromagnetic fields. This approach is more flexible than chem- 
ical and mechanical approaches. Not only electromagnetic fields are easy to control, but being tailored properly, 
they enable the creation of graphene-based tunable elements, including quantum wires and waveguides with unique 
properties, such as unidirectional conductivity, robustness to disorder, etc. The diversity of methods and approaches 
makes it increasingly difficult to summarize of the current state-of-the-art in this area, and calls for a systematic 
classification by both methods and results. An attempt of such classification is a goal of this chapter. 

Manipulating charge carriers by electric fields (i.e., addin g sca lar potentials of different shapes to the Dirac equation) 
is a very popular approach (see, e.g., the reviews in Refs. 0,1207}). This method has provided interesting and surprising 
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outcomes (i.e., |4l. l207H209l |). Unfortunately, it turned out to be impossible to repel or localize electrons in all directions 
and at all energies by only using an electric field [|[ . There is always a channel in any electric- field barrier where the 
charge can escape through. 

The charge-confining and current-guiding capabilities produced b y ma gnetic barriers are well known and have 
already opened certain possibilities for practical applications (e.g., (210l |). However, when it comes to designing 
fast-tunable electronic devices (switches, filters, etc.) a difficulty emerges: most of the existing magnetic-barrier 
technologies usually involve the deposition, either on top or beneath the graphene sheet, of a pattern of magnetic 
material, which reproduces the desired magnetic field distribution in the sample. Any subsequent change of parameters 
would require building a new setup, creating formidable (if surmountable) obstacles for harnessing magnetic barriers 
as elements of fast-acting electronic devices. Grathene structure s based o n the effective magnetic field created by 
applying inhomogeneous pressure or strain to the graphene sheet [204| - |206| might also be useful for applications. 

An efficient way around this problem is the simultaneous use of inhomogeneous magnetic and electric fields. The 
proper combination of these two not only preserves (or even improves) the necessary transport properties of graphene 
samples, but makes them easy to control by tuning the spatial distribution of the electric potential, for a fixed magnetic 
field. 

Hereafter, we focus on low-energy excitations when the inter-valley scattering |21l| is negligible and quasiparticles 
can be considered as massless Dirac fcrmions. To simplify the presentation, all spin-related effects are also neglected 
below. 

The building block of all charge-confining elements is a field- induced barrier (i.e., a reflecting wall). To introduce 
this in the most general way, consider two graphene half-planes, x < xq — 1/2 and x > xo + 1/2, subject to different 
stationary electric V* and magnetic A, potentials (i = 1,2), as shown in Fig. [23 Here we assume that A = A(x)y, 
which means that the magnetic field B — B(x)z. When \x — xq\ < 1/2, the fields and potentials are x-dependent, 
providing the transition from {Vi,A±,Bi} to {V2, A2, B2} (see Fig. [231) . If the width I is large, compared to the 
graphene lattice spacing a, and much smaller than the Fermi wavelength A^, a <C I <C Xp, then the smooth profile 
can be replaced in the calculations by a step function as shown by the dotted line in Fig. (23 When B = we have 
pn junction, which is discussed in Section VI. 



P 




FIG. 25: (Color online) Depending on the problem being considered in this section, the function P(x) is either one of the 
potentials V(x), A(x), or the magnetic field B(x). Here, P(x) is shown as a continuous green line, while the red dashed line 
repeats the step-like approximation. 



A. Magnetic barriers: V = 

As it was demonstrated theoretically in Refs. |l50l 12121 |213| . charge carriers in graphene could be confined in 
space by applying a properly tailored static inhomogeneous magnetic field perpendicular to the graphene plane (a;, y). 
There are numerous unusual and intriguing features in this phenomenon, but at its basic level, this happens because 
the trajectory of the quasiparticle incoming from free space x < bends inside the magnetic barrier (Fig. 125 j) and 
eventually exits backwar d, in dependently of its cyclotron radius, thus making the wall a perfect reflector for Dirac 
electrons of any energy |214| . Obviously, a barrier of finite width d reflects only charges with energies bellow a 
threshold, namely those whose cyclotron radius is smaller than d. Below we consider two types of these magnetic 
barriers. 
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1. Step-wise vector potential A: S -function-like magnetic field B 



This basic magnetic field profile, B(x) = Bq £b S(x) (here £b = eB§ is the magnetic length), already displays 
an important advantage of magnetic barriers over electric ones: there is a ran ge of energies for which it is a perfect 
mirror for Dirac electrons for all angles of incidence including perpendicular [210|. Outside this energy range, the 
transparency is angular-dependent, and the analog of the total internal reflection (TIR) of light takes place when 
the angle of incidence exceeds some critical angle </>tir, which depends on the applied magnetic field. The resulting 
charge refraction is not of t he c lassical Snell's law type, and is asymmetric with respect to the transformation of the 
angle of incidence — > — 4> |210j | . 

Rather counter- intuitive is also the fact that a "point-like" magnetic field, B(x) oc S(x), creates a bound (in the 
x-direction) state at zero energy. This state is a "linear" analog of two-dimensional surface waves: it is exponentially 
localized in the x-direction around the line x = 0, and propagates along the y-axis with wave numbers ± k y . The 
longitudinal wave number k y is enclosed in the finite interval £^/2), and the x-dependence of the wave 

function is asymmetric, so that the r ates of the exponential decay at x < and x > are different, depending on the 
sign of k y (direction of prop agati on) |214| . The possibilities to experimentally produce highly localized magnetic fields 
have been discussed in Ref. [210l ] , where charge transport in the presence of various arrangements of 5- type magnetic 
barriers was studied in terms of "electron optics" . 

2. Step-wise magnetic field B: piece-wise linear vector potential A 

The transport properties of Dirac electrons in a step-wise magnetic field [i.e., B(x) = Bi9{— x) + B 2 0(x)} are 
different and of greater variety than in the step- wise potential considered in the previous subsection I VII ATI While in 
the later case, the system is obviously invariant with respect to a shift of the A-stcp [namely, A(x) — > A(x) + const], 
for a step- wise B this is not true: if B\ and B 2 are parallel (7 = B\jB 2 > 0) or antiparallel (7 < 0) makes a significant 
difference. 

To better understand the behavior of Dirac electrons in nonhomogeneous magnetic fields, it is illuminating to 
compare it with the conventional two-dimensional electron gas. Fundamental differences already exist in the uniform 
field. In contrast with the classical Landau levels £ n in infinite space, £ n oc (n + 1/2), n > 0, the quantization of the 
cigen-energy of the Dirac equation produces £ n oc sign(n) \f\n\, with —00 < n < 00. Positive values of n, n > 0, are 
associated with electron-like charge carriers, while n < corresponds to holes. The cigenstates with n^0 are similar 
to the states of the conventional two-dimensional electron gas, whereas the zero-energy state n = possesses different 
properties. Bound states associated with this n — Landau level have different features than the states associated 
with n =/= 0. 

The energy in homogeneous magnetic fields does not depend on the wave number k, therefore the group velocity 
is zero, v g = d£ n /dk = 0, and the states carry no current. A dispersion, and therefore a non-zero group velocity 
v v = d£ n /dk y ^ can appear either for states localized near the boundary of a finite sample (so-called edge states), 
or due to the spatial inhomogencity of the magnetic field, X7 X B 7^ 0, which creates bound states localized in the 
x-dircction and propagating along th e ?/-a xis with its drift velocity proportional to V B x B (denoted, by analogy with 
edge states, as magnetic edge states [215j ). 

For the particular case Bi = (free half-space for x < 0, constant magnetic field B 2 > at x > 0), there is an 
infinite number of bound (in the x-direction) dispersive states labeled by the Landau-level index n, whose energies 
are proportional to sign(n) •y/jnj'. These states are localized as functions of x, centered around points whose locations 
depend on the wave number k y . Remarkably, these localized states exist only with one sign of the wave number k v , 
cither positive or negative, depending on the orientation of the magnetic field B 2 . Since the group velocity has opposite 
signs for +n and — n, the direction of the charge flow created by ele ctron s and holes is the same, and therefore any 
bound n-state carries a finite unidirectional current along the y-axis |210l | . Therefore, a particle in such a state never 
undergoes backscattering; hence it is practically insensitive to disorder and Anderson localization never takes place, 
no matter how strong the disorder. There is also a bound state with £ = 0, when k x = ik y . However, this state is 
dispersionless an d do es not carry any current. 

When Bi ^ |216l | and it is parallel to B 2 (7 > 0), the Landau levels of Dirac quasiparticles at large positive k y are 
localized in an effective potential well around x = ~k y £ B . With k y decreasing to negative values, the dimensionlcss 
energy levels gradually change to sign(n)y7jnj at large negative n and shift in space to x = —k y £ B /^. It is important 
to note that the directions of the drift (signs of d£ n /dk y ) are opposite for electron- and hole-like particles (±n), 
thus providing a non-zero total current. In the vicinity of k y = 0, the corresponding states become localized at 
x = 0. Similarly to B\ = case, the magnetic barrier with 7 > supports bound states, which create unidirectional 
conductivity at n 7^ 0, while zero-energy solutions carry no currents. 
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When 7 < (for antiparallcl B\ and B 2 ) [216j | . the effective potential at large positive k y has two minima (two 
connected harmonic wells), located far away from the boundary x = 0, so that the states are localized in each well. 
As k y moves to negative values, the effective potential shifts toward x = and transforms into a single non-harmonic 
potential well. The eigen-energies with n correspond to states which support non-zero unidirectional current 
following classical, so-called snake, orbits [2 1 51 ] confined to a na rrow one- dimensional channel centered at the line 
x = where the magnetic field changes its sign. It is shown in [2 1 7t 12181 ] that in a symmetric graphene sample of 
a finite width L (-L/2 < x < L/2) this current is compensated by real edge states localized close to the sample 
boundaries. The states with n = exhibit both electron and hole features, which is highly unusual and is unique for 
Dirac quasi-particles in graphene. 

B. Combined magneto-electric barriers: B 7^ 0, V ^ 

In principle, the charge-confining and guiding capabilities of magnetic walls presented above open up certain pos- 
sibilities for practical applications. However, as it was mentioned before, the parameters of the magnetic barriers 
cannot be changed fast enough, which makes it problematic to use them as elements of fast-acting electronic devices. 

To overcome this problem, it is convenient to use a combination of inhomogeneous magnetic and electric fields, 
which enables the efficient control of the transport properties of graphene samples by tuning the electric potential 
without changing the parameters of the magnetic field. 

To introduce a basic setup combining magnetic and electric fields, we now consider a single magneto-electric barrier 
produced by superimposing a scalar potential V of the same step-like shape on the magnetic structure with a (5-like 
magnetic field [i.e., a step-wise vector potential A(x) = A\6{— x)+ A 2 9{x)]. This system possesses unique properties 
that make it different from other types of barriers. Graphene subject to mutually perpendicular electric and magnetic 
fields supports states which are localized near the barrier. These curre nt-carrying states (surface waves) correspond 
to quasiparticles moving along the barrier only in one direction |219l ]. This direction, as well as the value of the 
quasiparticle velocity, are easily controlled by the electrostatic potential. These states correspond to the classical drift 
of charged particles in crossed electric and magnetic fields. They exist if and only if the drift velocity Vd = cE/B is 
smaller than the Fermi velocity vf [here E ~ (V 2 — V\)/l and B ~ (A\ — A 2 )/l are the electric and magnetic fields in 
a barrier of finite length I] . The ab sence of counter-propagating states prevents the backscattering induced by either 
irregularities in graphene [l88l l208| ] or by the fluctuations of the magnetic field. 

For potential applications, the important feature of a single magneto-electric barrier is that the transport (electric 
current) across or along this structure can be controlled by manipulating only the electric potentials V\ and V 2 - In 
particular: 

— The transmission and reflection coefficients across the junction between two areas with different values V\ , A\ 
and V 2 ,A 2 (Fig. [25]), and the angle of refraction (i.e., the direction of the transmitted current) depend on the electric 
potentials. Specifically, tuning V\ and/or V 2 can change the angle of incidence where the barrier is totally transparent, 
and thus the Klein tunneling can be observed; 

— When the inequality 

\V 1 -E\ + \V 2 -£\<\A l -A 2 \ (5) 

holds (here £ is the energy of the quasiparticles; all units are dimensionless), the step is a perfect reflector for electrons 
at all angles of incidence and the junction is locked for the electric current. 

— If 

\V 1 -V 2 \<\A 1 -A 2 \, (6) 

a wave (current) exists which propagates unidirectionally along the barrier (in the y-direction) with the dimensionless 
group velocity v = Vd/vF < 1 and is exponentially localized in the x-direction. 

C. Waveguide with electrically-tuned parameters 

While one barrier forms a wire, two such barriers constitute a waveguide. This waveguide supports modes that 
are similar to the electromagnetic eigenmodes of a dielectric waveguide and likewise have quantized transverse wave 
numbers (The ana logy between t he t ransport of Dirac electrons in graphene and light propagation in dielectrics is 
described in Refs. [1791 Il88l 1 1991 . |219| ). However, along with them there is another set of waves that is appropriate 
to call "extraordinary" [219(. They are formed by two coupled surface waves propagating along the waveguide walls 
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(barriers). There is an energy gap where only extraordinary modes exist. Decreasing the spacing between the barriers 
broadens this gap. The extraordinary modes are also stable against backscattering. 

An important feature of field-induced waveguides in graphcne, which is favorable for the creation of tunable elec- 
tronic devices is that the transport properties of these structures are strongly dependent on the parameters of the 
barriers. These parameters are the potentials Ai tT and Vi, r of the left and right semiplanes, respectively, surrounding 
the central region where the potentials are equal to zero, A c = V c = 0. In particular: 

— When 

Ai = A r and Vj = V r 

and the inequalities (|5|) and (|6|) are valid for both barriers (waveguide walls), then the extraordinary modes are 
unidirectional. This makes them immune to backscattering, and therefore robust against y-depcndcnt disorder. 
- If 

Ai = -Ar and Vi = V r , 

then the surface waves "attached" to the barriers propagate along the y-axis in opposite directions and the extraor- 
dinary modes are bidirectional. Nevertheless, the backscattering also does not affect the total current, due to the 
spatial separation of the charge fluxes with opposite directions. 

— When 

Ai = -A r and Vi = —V r , 

the spectrum of the extraordinary modes is independent of the distance between the barriers, and therefore there is 
no cutoff energy for them. This means that extraordinary modes can penetrate through an arbitrary narrow part of 
the waveguide. 

In Tables [X] and IXI1 the barriers and waveguides of the above mentioned types and combinations of fields are 
categorized according to the following features: 

• Their ability to reflect all incident current (perfect wall) 

• Their ability to support bounded electron states 

• The type of spectrum of the propagating modes (either continuous or discrete) 

• The directionality of the current (either uni- or in bidirectional) 

• A separate column lists some other distinctive features because graphcne in magnetic fields has unusual prop- 
erties. 

In summary, changing the electric potential (with the magnetic field unchanged) one can switch on and off the 
current through the barrier and create/destroy a unidirectional quantum wire along the barrier. Moreover, changing 
the electric potential one can create waveguides with unique, exotic transport properties. 
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Barrier type 


Per- 
fect 
wall 


Bound 

states 


Directionality of the 
bound state current 


Comments 




V 


Electric 
potential 

1 Vl 

X 


No 


No 


No bound state current. 


When min{K} < £ < max{K}, the 
optics analogy of the barrier is the in- 
terface between two dielectrics with op- 
posite signs of the refraction indexes. 
Otherwise, this is the interface between 
usual dielectric media. Total inter- 
nal refraction is possible. Cannot be 
opaque for all angles of incidence (i.e., 
Klein tunneling). 


A X .-.. 


A 


Vector 
potential 

Al 

X 


Yes 


Yes 


One bound state with zero 
energy {£ = 0) and zero 
group velocity (v g = 0) 
along the barrier, and, 
therefore, carries no cur- 
rent. This bound state is 
associated with the Lan- 
dau level with n = 0. 


In a certain range of energies the bar- 
rier is opaque for all angles of incidence. 
This barrier is similar to the barrier 
generated by a graphene sheet strain. 
The difference is that the strain gener- 
ates an effective vector potentials jump 
(effective magnetic fields) with oppo- 
site (due to the time-reversal symme- 
try) signs in two valleys, whereas the 
real magnetic field has the same sign in 
both valleys. 


B\ 


B 


Ivldyi IcLIU 

field 

B2 

X 


Yes 


Yes 


Bi-directional conductivity 
when both Bi and B2 
are non-zero; otherwise 
uni-directional conductiv- 
ity along the barrier. 


When 7 = B1/B2 > 0, the bound 
state is similar to the classical elec- 
trodynamics state with X7B x B drift. 
When 7 < and Landau-level index 
n / 0, the state is similar to a snake 
state (charged particle motion along 
the B = line). Can be opaque for 
all angles of incidence. 


Ay... 
V\ — 


V,A 
potei 


1 IO dl IU VoL/LUI 

itials 

Vl 

A 2 X 


Yes 


Yes 


Confined uni-directional 
state with linear spec- 
trum d£/dky = Vd when 
\v d \ = c|(V a - Vi)/(Ai - 
A2)| < «F- In a certain 
range of energies, the 

Ucll 1 lei la UJJcLL^Llc 1U1 aLL 

angles of incidence. Easily 
controlled by the electric 
potential. 


The classical electrodynamics analogy 
of the bound state is the charged par- 
ticle drift in crossed electric and mag- 
netic fields. The drift velocity is Vd = 
cE/B = c{V 2 - Vi)/(Ai - A 2 ) 


v\— 


V, B 
biec 
and 


trie potential 
magnetic field 

v 2 

B 2 X 


Yes 


Yes 


Same as the magnetic field 
barrier shown above. 


Same as the magnetic field barrier 
shown above. 



TABLE X: Basic properties of electric and/or magnetic barriers in graphene. The ability to reflect all incident electric currents 
(or the ability to be a perfect reflecting wall) is indicated in the second column. The presence or absence of bound states is 
listed in the third column. Additional physical properties are summarized in the remaining columns. 



35 



Waveguide type 



Guided modes 



Comments 



Discrete set of bulk waves. 



Optics analogy: this is a dielectric 
waveguide based on total internal 
reflection. 



Vy 



Discrete set of bulk waves 
and one surface mode with 
zero group velocity along the 
waveguide. 



- A r 



The bulk waves are similar to di- 
electric waveguide modes. The 
surface mode is localized near the 
waveguide walls and consists of two 
coupled surface modes associated 
with the barriers (the waveguide 
walls). 



5/ 



B r 



Discrete set of bulk waves 
when a = Bi/B r = ±1 and 
B c = 0. The mode is lo- 
calized near the boundaries 
when a — 1 and B c < 0. 



The discrete set of bulk waves 
is similar to the eigenmodes of 
a dielectric slab surrounded by a 
medium with negative permittiv- 
ity/permeability. It always has a 
snake state with noncompensated 
current. When a = 1 and B c < 0, 
there are counter-propagating cur- 
rents near the walls (two snake 
states near the walls). 



Vl 



V, A 



Discrete set of bulk waves 
and two surface modes. 



Bulk waves are similar to the eigen- 
modes of dielectric waveguides. 
When the spacing between the bar- 
riers is rather small, there are only 
surface modes. The propagation 
direction of surface modes is con- 
trolled by the electric potentials. 



x 

- A 



r 
V r 



TABLE XI: Basic properties of waveguides in graphene. These are produced by various configurations of the applied electric 
and magnetic fields. Note that the convention used here for dotted and dashed lines is different from the one used in the 
previous Table. 



36 



VIII. CONCLUSIONS 

Graphene is a material with many interesting features which make it an attractive candidate for microelectronic 
and micromcchanical applications. In the above pages we very briefly outlined several ideas and notions driving 
current graphene mcsoscopic research: edge states, geometric quantization, quasi-bound states, Coulomb blockade, 
Klein tunneling through pn-junctions, etc. Some of them, like Klein tunneling, are unique to graphene. Others, e.g., 
Coulomb blockade and geometric quantization, have a much longer history. Yet, even in the latter case, the peculiar 
properties of graphene give rise to new features, for example, a much larger energy scale for the confinement inside a 
quantum dot. 

Although, many of the theoretical studies of graphene mesoscopic systems are done in the single-electron approxima- 
tion, the use of many-body techniques are often warranted. Indeed, the single-electron approximation could introduce 
qualitative errors, which may be corrected only if proper many-body effects are accounted. 

Numerous theoretical proposals have not been explored experimentally. For instance, the realization of nanoribbons 
with atomically-sharp edges remains a distant possibility. Many suggested devices impose stringent conditions on 
samples, in terms of purity and regularity of the sample geometry. Some of these proposals are stimulating various 
current experiments. An important direction in this quest is to control and to understand the disorder, which may 
enter through many routes: as foreign atoms adsorbed on samples or chemically attached to the edges, as imperfections 
of the sample edges, as random elastic deformations, or as bulk defects of varied nature. On the other hand, disorder 
is not always an enemy, since its use may be beneficial under certain circumstances. 

Graphene studies are still in their infancy, and it is too early to guess which of its unusual features will be more 
useful for applications. Yet, fabrication of several prototypic microelectronic devices, like field-effect transistors, 
biosensor, and integrated circuit, have been reported. In addition, graphene presents an excellent playground for 
fundamental condensed matter research, exciting enthusiasm of both experimentalists and theorists in numerous 
subfields of condensed matter physics. 
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Appendix A: Tight-binding model of graphene lattice 



In this Appendix we present basic notions necessary for the theoretical description of the single-electron properties 
of a graphene sheet. It provides some details omitted in Sec. [TTJ 

It is common to describe a graphene sample in terms of a tight-binding model on the honeycomb lattice. Such 
lattice can be split into two sublattices, denoted by A and B. The Hamiltonian of an electron hopping on the graphene 
sheet is given by Eq. ([!]), in which R runs over sublattice A 



m = 81 + aini + a 2 n 2 , 
where the primitive vectors of the honeycomb lattice are 

ai = ao(3/2,V3/2), 
a 2 = oo(3/2, -V3/2), 
and ni 2 arc integers. The vectors Si (i = 1,2,3) connect the nearest neighbors. They are 

Si = a o (-l,0), 

8 2 = a (l/2,V3/2), 

S 3 = 00(1/2,-73/2). 

The geometry of the graphene lattice is presented in Fig. [1] 
The Schrodinger equation can be written as 



i=l,2 



(Al) 

(A2) 
(A3) 

(A4) 
(A5) 
(A6) 



(A7) 
(A8) 



i=l,2 



where tp£ (tpR+sJ denotes the wave function value at the site R (at the site R + Si) of sublattice A (sublattice B). 

The primitive cell of graphene contains two atoms, one at R, another at R + Si. Therefore, it is convenient to 
define the two-component (spinor) wave function 



R 



4> 



By construction, the function is defined on sublattice A, Eq. (|A1[) . 
The action of H on a plane wave 



*r = *k exp(-ik • R) 



can be expressed as 

tk = t 

For every k there are two eigenstates 



1 + 2exp 



.3k x a 



-tt 







kydQ 



^k± 



with eigenvalues: 



exp (i6» k ) 



e k± = ±\t k \=±ty/3 + F(k), 

'V3 



l*k| 



F(k) = 4 cos ( -k x a Q ) cos 



2 cos 



(A9) 

(A10) 

(All) 
(A12) 

(A13) 
(A14) 

(A15) 
(A16) 



38 



The states with negative (positive) energy are filled (empty) at T = 0. The allowed values of k lie within the Brillouin 
zone presented in Fig. [2] 

The reciprocal lattice is characterized by the following lattice vectors 

d! = (47r/3a ,0), (A17) 
d 2 = (-27r/3a ,27r/V3a ). (A18) 

The amplitude tk and energy £k± are invariant under shifts over di i2 . 

The quantity £k± vanishes at the six corners of the Brillouin zone: K h2 = (0,±47r/(3\/3ao)) and K 3 , 4:5i6 = 
(±27r/(3a ),±27r/(3V3a )). These are the locations of the famous Dirac cones of graphcne. These six cones can be 
split into two equivalence classes: all cones inside a given equivalence class are connected by a reciprocal vectors. We 
choose cones located at Ki j2 to be representatives of these classes. 

If one is interested in the low-energy description only, then the tight-binding Hamiltonian may be replaced by the 
Dirac Hamiltonian with the help of the following derivation. Near the cones, the Hamiltonian Eq. (|A11|) may be 
expanded in orders of k — Ki j2 . The spinor itself may be represented as 

*(R) = *i(R) exp(-iKi • R) + * 2 (R) exp(-iK 2 ■ R), (A19) 

and the Schrodinger equation splits into two copies of the Weyl-Dirac equation ([3]), where we treat 5 , i 2 (R) as slowly 
varying functions of the continuous variable R, and the plus (minus) sign is chosen for "Ji (vl' 2 ). 

Describing the electronic properties of the graphene in terms of Eq. ([3]) is very popular for two reasons. First, it 
is much simpler than the full tight-binding Hamiltonian ([IJ. Second, unlike the tight-binding Hamiltonian, which 
may require knowledge of several hopping amplitudes, only one parameter v-p has to be specified. Of course, the 
tight-binding model is needed when a more detailed description is required. 



Appendix B: Edge states 

Usually, two types of edges are discussed in the literature: zigzag and armchair. A variant of the zigzag is the Klein 
edge [H|]. All three types are shown in Fig.Q] 



1. Armchair edge 

The physics of the armchair edge is simpler than that of the zigzag edge because the latter supports zero-energy 
localized states, while the former does not. The easiest way to describe an electron near an armchair edge is to use the 
Weyl-Dirac equation ([3]) with the appropriate boundary condition. The general problem of the boundary condition 
for the Weyl-Dirac equation is investi gate d in Rcfs. (3~H. [36l l50l - l52j . Here we use a simple explicit form of the boundary 
condition suitable for armchair edge [811 ] . Namely, we demand that our spinor wave function Eq. (|A19[) vanishes at 
the edge, which we assume to be located at y = 

*i(R)|„=o = -*2(R)|„=o. (Bl) 
For an infinite half-plane, the solution of the Weyl-Dirac equations with this boundary condition is equal to 

*i±(R) = -7= [ \ h wj.) exp(-ik x x - ik y y), (B2) 



2 \T(*^ + k y )/k 

* 2±(R) = ~7! (t(^+*v)/0 ex PH^ + ^)- (B3) 

The eigenfunctions \&i,2+ (^1.2-) correspond to positive (negative) eigenvalues. The total wave function is to be 
constructed according to Eq. (|A19|I . 

However, it is not necessary to use the approximate description in terms of the Weyl-Dirac equation. The tight- 
binding model may be solved near the armchair edge as well. To construct such a solution, note first that the atoms at 
the very edge are special: they have only two nearest neighbors, while the atoms in the bulk have three. This means 
that Eqs. (|A7|A8j) must be modified to account for this fact. It is more convenient, however, to introduce auxiliary 
rows of carbon atoms (Fig. I26|) and demand that the wave function vanishes on these additional atoms. Since there is 
no wave function density on the auxiliary atoms, they do not contribute to the Schrodinger equation for the physical 



.V 
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FIG. 26: Introduction of auxiliary atoms (hatched circles) at an armchair edge. 



atoms at the edge. This construction is used in several papers (e.g., Ref. [8S 
is 



*r = *k± exp(-ik x x) sin[k y (y + V3a /2)}. 



1901). The solution of Eqs. £57) and 



(B4) 



This wave function vanishes at y = — v3<Zo/2, which is where the auxiliary atoms are located. At the physical edge 
of the nanoribbon (y = 0), however, the electron density remains non-zero. Note that momentum components k x>y 
in Eq. (|B4[) is measured from the center of the Brillouin zone. However, in Eq. (|B2[) and Eq. (|B3[) they are measured 
from the Dirac cones. When this subtlety is accounted it is easy to demonstrate that both solutions have the identical 
dependence on R. They differ only in the value of the spinor part ^k- This is because Eq. (|B2p and Eq. (|B3[) are 
only approximations which are accurate at small energy only. 

The discussion presented above assumes that the properties of the carbon-carbon bonds near the edge remain the 
same. In a real sample t his a ssumption is only an approximation. A variety of perturbations may be present near the 
edge. For example, Ref. [220j ] investigates the edge structure with the help of density-functional methods. It is found 
that the carbon-carbon bonds at armchair edges are shorter than in the bulk (1.26 A versus 1.4 A). This means that 
the effective electron hopping at the edge t c differs from its value t in the bulk. 

Also, non-carbon radicals may be attached to the unpaired chemical bond of carbon atoms at the armchair edge. 
This means that, depending on the description, either the Hamiltonian near the edge, or the boundary conditions 
have to be modified. Fortunately, these alterations are small, and may be treated as weak corrections. 



2. Zigzag edge 

The physics of the zigzag edge is richer than the physics of the armchair edge due to the zero-energy states localized 
at the edge. When graphene is described by the Hamiltonian (p}, which only contains nearest- neighbor hopping, these 
states are dispersionless and macroscopically degenerate. 

The simplest way to detect the presence of edge states is to use the Weyl-Dirac equation with the boundary condition 
appropriate for the zigzag edge [HI, [H, H(| ■ For t he p urpose of demonstrating the presence of edge states the simplest 
version of the boundary condition is used here [8 if . The zigzag edge, unlike the armchair edge, consists of atoms 
belonging to the same sublattice. This is clearly seen in Fig.[TJ all atoms at the left edge belong to A sublattice (red). 
The right zigzag edge has all its atoms on the B sublattice. Extending the left edge by a column of the auxiliary 
atoms, we demand that the wave function vanishes on them (assume that the auxiliary atoms are located at x = 0): 



?A B U=o = 0. 

Then the following spinor functions are the solutions of the Weyl-Dirac equation: when k y > 



(B5) 



*i(R) 



exp(—k y x + ik y y), * 2 



(B6) 
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and when k y < 



*i = 0, * 2 (R)= (ljexp(-\k y \x + ik y y). (B7) 

Both solutions decay for large values of x and correspond to the zero eigenvalue. The difference between them is that 
Eq. (|B6[) describes the solution near K l5 while Eq. (|B7|) near K 2 . The wave function Eq. (|A19[) constructed from 
Eqs. (|B6| and (JBTJ) is equal to 

*(R) = (^j cxp(-iKi • R - k y x + ik y y) + ( * j cxp(-iK 2 ■ R — k' y x — ik' y y), (B8) 

where both k y and k' y are positive and small compared to 1/clq. The quantity k y sets a parameter with dimension of 
length 

Acdgc = | (B9) 

which characterizes how deeply the edge states extend into the bulk of graphene. 

Since the Weyl-Dirac equation is only an approximation of the tight-binding Hamiltonian (JlJ, we cannot reliably 
use Eq. ()B8[) for large values of k y and k' y . However, the tight-binding problem may be solved exactly [HI to discover 
that there is a degenerate manifold of states labeled by the momentum k y , which stretches from Ki to K 2 . 



a. Effect of the longer-range hopping 

The degeneracy of edge states is purely accidental property of Hamiltonian ([1]). It disappears when other terms are 
added to H. For example, in Refs. [5514581 the effect of next-to-nearest-neighbor hopping on edge states is studied both 
analytically and numerically. In Ref. [581 ] the authors made three keen observations: (i) the next-to-nearest-neighbor 
hopping effectively induces a shift in the local potential; (ii) the shift depends on site's position, namely, near the 
edge it is not the same as in the bulk; and (Hi) this spacial variation of the potential induces finite dispersion for an 
otherwise dispersionless edge states. To prove that (Hi) holds true, consider the following argument. The states with 
large A e dge are insensitive to the potential variation (they "feel" the potential averaged over large A ec ige), however, 
those with small A e d ge are affected strongly; since A e( jge depends on \k y \, the edge states acquire the dispersion. 

This discussion suggests that, if we were to describe this phenomena with the help of Weyl-Dirac equation, the 
boundary conditions, Eq. (|B5[) . must be modified to account for the effect of the potential modulation near the 
boundary. Indeed, it is demonstrated in Ref. [52j that one can generalize Eq. (|B5|) and reproduce the findings of 
Ref. (58|, at least near the apexes of the Dirac cones. 
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